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SUMMARY 


The nu) s L noticeable effect of air pollution «m t lie ptopert Iva 
of the atmosphere is the reduction in visibility. J3y consult ting 
dynamical and microphysical processes, the present study invest i gates 
the formation of adveetion fog associated with aerosols, due to 
combustion-related pollutants, as condensation nuclei. We are 
particularly interested in the study of the effects of mult l-dispersi 
distribution of the aerosol population, and also single component and 
multi-component aerosol species on the condensat ion/uuc Leat ioti e.ocesse:, 
which affect the reduction in visibility. 

The results show that the aerosol population with a high particle 
concentration provides more favorable conditions for the formation oi 
a denser fog than the aerosol population with a greater particle sice 
distribution if the value of the mass concentration of the aerosols is 
kept constant. 

The results of the present study can be used as numeric ul predict U-n 
of fog formation. Two-dimensional observations in horizontal and vert ical. 
coordinates, together with time-dependent measurements are needed as 
initial values for the following physical parameters: (1) wind prof iles? 

(2) temperature profiles; (3) humidity profiles; (4) mass concentration 
of aerosol particles; (5) particle size distribution of aerosols; and (h) 
chemical composition of aerosu.s. Formation and dissipation of adve.u um 
fog, thus, can be forecasted numerically by introducing initial values 
obtained from the observations. 

The present research supported G. S. Liav, a former graduate 
student, who obtained his degree of Doctor of Philosophy from The Uaiv.a.i 
of Alabama in Huntsville. 
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Figure 1. Three typos of aerosol population diotribut ions. 

Kach distribution is a combination of Lin cv step 
function;; in which the total mans; eoncentrat ion 

conserved and equal to 54. V ug/m^. 1H 

Figure 2. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, ,,SO, , 

with type A aerosol population distribution, an 
shown in Figure i, at two horizontal locations, 
x - 3 and 6 km, and both at vertical location, 
z ~ lm. 21 


Figure 3. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, (NH^).,SO^, 

with type A aerosol population distribution, as 
shown in Figure 1, at two horizontal locations, 
x = J and 6 km, and both at vertical location, 
z = 4 m. 

Figure 4. Visibility and relative humidity profiles of fog 

formation due to condensation nuclei, (Nil, ) J1U, , 

4 2 4 

with type B aerosol population distribution, as 
shown in Figure 1, at two horizontal locations, 
x = 3 and 6 km, and both at vertical location, 
z = 1 m. 

Figure 5. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, (Nil,),, SO,, 

with type B aerosol population distribution, as 
shown in Figure 1, at two horizontal locations, 
x ~ 3 and 6 km, and both at vertical location, 
z = 4 m. 


Figure 6. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, (NII,),,Su,, 

4 2 4 

with type C aerosol population distribution, as- 
shown in Figure 1, at two horizontal locations, 
x - 3 and 6 km, and both at vertical location, 
z = 1 m. 

figure 7. Visibility and relative humidity profiles of fog 

formation due to condensation nuclei, (NH,)_SO,, 

4 2 4 

with type C aerosol population distribution, as 
shown in Figure 1, at two horizontal locations, 
x = 3 and 6 km, and both at vertical location, 
z = 4 m. 
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Figure y. Aerosol population distribution wiLh a fixed value 
of maos concentration: (a) High particie concen- 
tration - shifted dial ribuLitm; (b) Low particle 
concentration - shifted distribution. 

Figure 9, Three types of aerosol population distributions. 

Kaeh distribution is a coubin.it ion of three step 
functions in which the total mass concentration 
are conserved and equal to 54.9 ;.g/ m^, for 
(NH^).,S0 4 , and iy.3hg/m3 for Ca(NO J ) 2 . The solid 

line denotes the distribution for and 

the broken line implies the distribution for 
Ca(h'0 3 ) 2 . 

Figure 10. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, (NH, ) 2 S0^ 

and Ca(N0-j) 2> with type A aerosol population 

distribution, as shown in Figure 9, at two 
horizontal locations, *; = 3 and o km, and both 
at vertical location, a •- lm. 

Figure 11. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, 

and Ca(N0 3 ) 2 , with type A aerosol population 

distribution, as shown in Figure 9, at two 
horizontal locations, x ~ 3 and 6 km, and both 
at vertical location, z ~ 4m. 

Figure 12. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, (NH^) 2 30^ 

and Ca(N0 3 ) 2 , with type B aerosol population 

distribution, as shown in Figure 9 , at two 
horizontal locations, x - 3 and 6 km, and both at 
vertical location, z = lm. 

Figure 13. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, (NH^) 2 S0^ 

and Ca(N0 3 ) 2> with type B aerosol population 

distribution, as shown in Figure 9, at two hori- 
zontal locations, x = 3 and 6 km, and both at 
vertical location, z = 4m. 

Figure 14. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, (KU^) 2 S0^ 

and Ca(N0 3 ) 2 , with type C aerosol population 

distribution, as shown in Figure 9, at two 
horizontal locations, x ~ 3 and 6 km, and both 
at vertical location, z = lm. 
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Figure 1J>. Visibility ..tni relative humidity profile-, of tog 
formation due to condens.it ion nuclei, 

4 it 

and Cu(N')..) , , with typo C aerosol population 

distribul ion, a:; , shown in Figure 9, at two 
horizontal location;;, x <= 3 and 6 km, and both 
at vertical location, z a 4m. 

Figure lb. Time rate of change of the amount of liquid water content 

associated with condensation nuclei, (N'll.J.dh), and I'aCN’u j,., 

4 2 a 5 l 

with type A aerosol population distrihut ion, as shown in 
Figure 9, at two horizontal locations, ~ 3 and b km, and 
both at vertical location, s = 1 m. 

Figure 17. Time rata* of change of the amount oi liquid water content 

associated with condensation nuclei, (NH. ) 0 S»), and OafNo ) ,, 

with type A aerosol population distribution, as shown in 
Figure 9, at two horizontal locations, x “ 3 and o km, and 
both at vertical location, z ~ 1 m. 


Figure 18. Time rate of change of the amount of liquid water content 

associated with condensation nuclei, (NH^),/!*)^ and daf.JJO ,) ,, 

with type 15 aerosol population distribution, as shewn in 
Figure 9, at two horizontal locations, x - 3 and b km, and 
both at vertical location, z = 1 m. 

Figure 19. Time rate of change of the amount of liquid water content 

associated with condensation nuclei, (NIl^),,SO^ and CaCNP.^, 

with type 15 aerosol population distribution, as shown in 
Figure 9, at two horizontal locations, x - 3 and b km, and 
both at vertical location, z G 4 m. 


Figure 20. Time rate of change of the amount of liquid water content 

associated with condensation nuclei, and CaiAO.^ 1 a * 

with Type 0 aerosol population distribution, as shown in 
Figure 9, at two horizontal locations, x ~ 3 and b Ian, and 
both at vertical location z 1 m. 


Figure* 21. Time rate of change of the amount of liquid water content 

associated with condensation nuclei, (NH„ ) , ; SO /} and Ca(NO .),,, 

with type C aerosol population distribution, <,s .diown in 
Figure 9, at two horizontal locations, x - 3 and b km, and 
both at vertical location, z = 4 m. 
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I. INTRODUCTION 


Thu basic ingredients of a fop, <»r clued droplet a?'** 1 iquhl w at * i 
(or ice) and eondens.it Ion nuclei (or ice nuclei in the cast* of an ice 
fog or ice cloud). These nuclei are aerosol particles., e.eneraliv x.itl* 
diameters equal to or greater than 0.1 pm, and are nearly always at 
least partially hygroscopic. 

During the absence of aerosol particles in the atmosphere, con- 
densation of water vapor will occur only if the supersaturation roaches 
several hundred percent. However, in the atmosphere, there are abundant, 
foreign particles to serve as ready-prepared nuclei, both naturally 
produced and produced by energy-related combustion, for condensation 
and to present the large supersaturations from being achieved. 

The atmospheric aerosols cover a wide range of particle, si/e, 

-3 

from about 10 pm radius for the small ions consisting of a tew 

neutral air molecules clustered around a charged molecule to more lien 

10 pm for the largest salt and dust particles. Their concentration also 

3 

covers an enormous range, from less than 1000 per cm over the ocean, t. 
6 3 

perhaps 10 per cm in the highly polluted air of large industrial 

cities (Mason, 1371). A recent study of New York City aerosol particles 

4 

indicated that the average number of particles is (4 to 3) x 10 
-3 

particles cm for the lowest point at 5 AM, and for the highest point 
at 3 AM local time, with a sampling ranging in size from 0.1 to 1.3 |;m 
diameter (Leader er et al . , 1378). 

As to the chemical composition of atmospheric aerosols, the large 
quantities of SO^, NO^, MH^ and a strong acid such as II, ; S0^ have beet 
detected in highly industrialized areas (Leaderer et al . , 1978 ). A 


substant i ll portion ol thi:; aerosol product ion it; the result ol itiergv^ 
related i ik* 1 combust ion. Sulfur, an impurity iound both in coal ami 
petroleum, reacts with 0., at high temperatures producing ox idea ol S, 

80 . Furthermore, tiini’e air (79,’ N.,) it; the t)., t;ouree lor combust ion, 
the 1 orraat ion and subsequent emission ol NO , especially .m elevated 
combust ion temperature:;, are eontmon occurrences. The internal coin- 
bunt ion engine it; the dominant source ol the production ol NO,. This; 
production of SO., and NO due to energy-related Fuel combust ion in, 
in turn, transformed into aerosol particle.*! through photochemical 
reaction (Council on Environmental Quality, 1973). 

Based on the estimation made by the Council on Environmental 
Quality, energy conversion processes account tor 84"' of the SO , 97/, 
of the NO^, and 78% of the CO emissions (Council on Knviionmental 
Quality, 1973). However, the figures are quite different based on 
anthropogenic source points of view. For example, volcanic emissions 
account for more than 50% of SO.,, thunderstorms for mort m 75," 
of NO , etc. Since the dispersal of air pollution is unaffected by 
national boundaries, the emission of pollutants often results in 
global problems and concerns. 

Some of the atmospheric aerosols, in particular the salt aerosols 
with hygroscopic properties, are {mown to exhibit the property of 
deliquescence when exposed to a humid environment. The critical humidity 
at which a given salt particle abruptly transforms from the crystalline 
state to a saline droplet is related to the water vapor pressure in 
equilibrium with the saturated salt solution at the ambient temperature. 
Recently, Tang ej^ al. (1977, 1978) and Tang and Munkelwitx (1977), 
conducted a series of experiments to investigate the onset of 
growth which took place at a specific deliquescence humidity determined 
by the water activity at the eutonic composition lor some mixed salt 


aerosols NaGl-KGl, and (Nil,),, 80,-11,80, . The result showed that tin* 

onset of growth occurs at 73.8 + 0.32 relative humidity fur mixt d NaGl-KCl 
aerosols, and at 69.0 t 0.32 rolativi* humidity for sulfate aerot-xiis 
containing 0.75 to 0.95 molt* fraction (NH^SO^. 

In general, for a nucleus, or an aerosol particle to grow into a 
droplet, the air must have attained a certain degree of supersaturation 
(Fletcher, 1966; Weinstein and Silverman, 1973). At low supersaturation, 
only the nuclei large and hygroscopic, enough will go across the hump of 
the so-called Kohler curve (Mason, 1971). Once over the hump (i.e., 
the so-called critical size or critical supersaturation), these nuclei 
formally become, cloud droplets and can, in theory, grow to any size, 
provided a small amount of supersaturation is maintained. The figure 
for the presence of supersaturation to be a necessary condition to 
produce cloud droplets can be changed for condensation occurring in 
the polluted atmosphere (Hung and Liaw, 1979). 

The correlation of the fog formation with a high concentration of 
air pollutants associated with fuel combustion has been observed experi- 
mentally. With high concentrations of SO due to fuel combustion, the 

X 

formation of sulfuric acid by oxidation of SO^ on the surface of particles 
in a humid environment leads to the formation of small droplets of fog 
which would not otherwise have formed. For example, in determining the 
correlation between SO concentrations and fog formation in the German 
cities of Gelsenkirchen and Hamburg, Georgii (1969) found that 80 percent 
of all eases of high SO^ concentrations occurred with visibilities 
below 5 km. In addition to fog formation, increased precipitation has 
been linked to areas with high aerosol concentrations (Changnon, 1968). 


3 


The average mays concentration of aerosol particles from the Nat uni.il 
Air Sampling Network in the United States ranges from about 20 pg/m * in 
clean air in remote areas, to values of 100 to 280 pg/m * in urban areas 
(Hidy and Brock, 1971) . In the heavily polluted atmosphere of the hlgiilv 
industrialized areas, mass concentration with values as high as 2000 pg/m*, 

have been reported. As for particle number concentration, it varies from 

2 b J 

10 part iclas/<V: - ,r very clean air to 10 particles/cin in polluted air. 

The radius of ; an isation nucleus with combustion as a source ranges iium 

0.01 to 1 pm (Seinfeld, 1975). Thus, combustion is listed as a source of 

particles for sizes up to 1 pm diameter (Seinfeld, 1975) . 

The percentage distribution of mass concentration for SO is about 

18%, and for NU^ it is about 6% (Hidy and Brock, 1971) . This means that 

the mass concentration for SO is roughly less than 10 pg/nf* in clean air, 

and ranges from about 20 to 50 ug/in'^ in polluted air. The range ol mass 

concentration for NO is roughly one-’ third of the amount of the So range. 

X x 

Results obtained from Hung and Liaw (1980) show that the major 

contribution of combustion-related pollutants as condensation nuclei in 

the formation of fog comes from aerosols with species of SO . However, 

x 

N0 x as condensation nuclei also gives a considerable contribution to the 

formation of fog, but the weight is much less than that due to SO . In 

x 

this study, our major attention is focused on the effect of fog formation 
due to SO as condensation nuclei with the consideration of the influences 
made by the variety of particle size distribution. 

To elucidate the fog optical properties, Low et al. (1978) examined 
the effects of different spectral shapes an 1 widths on gamma and lognormal 
distributions with the same liquid water content, and concluded that the 
drop size distribution greatly affects the result of visibility even- 
though the liquid water content is the same. Observations show that 


4 


particle size distribution for combustion related aerosols are widely spread. 
In the present paper, we are interested in studying the effect of particle 
size distribution of aerosols as condensation nuclei for fog formation 
with keeping mass concentration fixed. To simplify the calculation, and 
also to make the results easy to compare, a combination of three step 
functions are chosen to modify the particle size distribution of aerosols. 

Acid rain is a major environmental problem on both sides of the 
Atlantic Ocean. Originally noticed and studied in the Scandinavian countries 
and in Canada, acid rain has been documented in this country; first in 
the Northeast and now throughout much of the United States. Increasing 
levels of acidity have already caused measurable damage to the environ- 
ment. Many lakes are now totally devoid of fish population; monuments and 
other man-made structures are being degraded; and yields of agricultural 
crops and forests may be decreasing. 

Acid rain, based on the definition of the terminology given by the 
Environmental Protection Agency, includes the wet removal processes such 

as rain and snow; the dry deposition component sucli as particulate sedi- 
mentation and/or gaseous aborption and adsorption; and special events 
such as fog, frost and dew which also contribute to the overall acidic 
deposition in amounts that need to be quantified. 

In the production of NO^, strong evidence indicates that the major 
producer comes from the automobile. There is no doubt that stationary 
source fuel combustions also contribute a heavy fraction of the total NO 
emitted t.o the atmosphere. Fuel-bound nitrogen compounds are ammonia, 
pyridine, and many other amines (Glassman, 1977). Similar to NO production 
by combustion, S as an inherent impurity in both coal and petroleum, SO 
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is also produced through fuel combustion. This production of 80 and N't),, 
duo to energy-related fuol combustion is, in turn, transformed into 
aerosol particles through photochemical reaction (Council on Environment al 
Quality, 1973). 

Energy conversion processes account for 842 of the SO , 97/ of the NO x , 
and 78% of the CO emissions eased on the estimation made by the Council on 
Environmental Quality (1973). However, the figures are quite different 
for data chosen from anthropogenic source points of view. I’or example, 
volcanic eruption, such as ash coughed up by Mount St. Helens, Washington, 
on May 18, 1980, could account for more than 50% of SU^, thunderstorms for 
more than 75% of NO , etc. Since the dispersal of air pollution is un- 
affected by national boundaries, the emission of pollutants often results 
in global problems and concerns. 

The production of SO and NO became the major sources of acid rain 
through the transformation of gas molecules to aerosol particles via photo- 
chemical reaction, and condensation processes in which the aerosol particles 
serve as condensation nuclei. The present study, in particular, deals 
with the mathematical simulation of the formation of advection fog and 
the accumulation of liquid water content in the atmosphere through the 
condensation processes in which the effects of the multi-disperse distri- 
bution of the initial aerosol population and multi-component species of 
aerosol particles as condensation nuclei are taken into consideration. 
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II. MATHEMATICAL MODEL 


A theoretical model is employed which describes the » volution of 
potential temperature, water vapor content, liquid water content., and 
horizontal and vertical winds as determined by the. processus oi vertical 
turbulent transfer and horizontal advection of momentum, energy and 
moisture, as well as radiation cooling, growth of water droplet;; based 
on microphysical processes, and drop sedimentation. 

The mathematical model is time-dependent two-dimensional in the 
X-X plane, and the boundary layer model, is assumed. All the quantities 
are uniform in the Y direction. The fundamental equations governing the 
macrophysical processes of the evolution of wind components, waLer vapor 
content, liquid water content and potential. temperaLure. under Lite in- 
fluences of vertical turbulent diffusion transfer, turbulent mu' alum 
transfer and turbulent energy transfer cm be expressed in the fullov.’iu;, 
three sets of conservation equations: 

(a) Mass conservation equation of the air and diffusion equations 
for water vapor and liquid water contents 
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liquation (2-1) is the expression of incompressible fluid for Lhe con- 
tinuity equation of dry air. This expression is Lrue only when Lhe wind 
velocity is much less than the speed of sound, in this case, Lhe wind 
velocity up to 100 m altitude is in general less Lhan 10 m s ^ which is 
less than 3% of the speed of sound. 
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The change of Lite water droplet radius, in Equation (.*'9) , due to 


t oudensation/ evaporation It; governed by mierophysieal processes, and 
given by (Cars tens et, al . , 1974) 


where 
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The nucluation process as considered in this paper, more specifically 
the constants of Equations (2-13) and (2-14), is valid only for pure 


hygroscopic salts. The overwhelming fraction or urban particulates is 
hydrocarbon coated or mixed and therefore much less active as a mid eating 
agent. Therefore, the proper modification is necessary for the model to 
be applied to the simulation of particular problems associated with urban 
pollutants as condensation nuclei for fog formation. 

In the present study, the diffusivity of water vapor in the air 
follows the model presented by Fuller et al . (1966) 


D = 


0.10! (Cf 1 + m' 1 f T 1 • 75 
~7ivy3 + v !/3 } 2 


1 1 . 


V 
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Tiu* other constants, employed In Equations (2-10) to (2-lh) which 
are temperature dependent cun he* simplified as iolluws (I'aluch, 19/lj: 


K = U. 1075 xUf h (T-273. 16) % 0.3/2a xlif' (eal/m= u K’/.ee; 


I. » 737.44 - U.52 T(oal./g) 


•J *75.7 - 0.148 (T ~ 273.16) (dynes era ), 


IM"; 

(.'-.Mi 


The semi-empirical expression is used for i t (*) and b, i.i ., 


P. , „ l v ) = 4,c 47 x 10 

l. <3 l 


■o ; 


.at 

5 3.;373 


• exp 


6/18, 


' 3 • it) T 


f 

■J 


k * p s. 3 t‘‘ u/ ' — 5.73o) .j 


Equation (2-20) shows how surface tension of the pure water depends 
on the temperature. A recent study of water drop.le.Ln exposed to polluted 
air indicates that the evaporation is greatly different over a surlaee 
with pure water than over a surface covered with a monolayer (Holier 
and Mallen, 1970). A Sky lab experiment involving water droplets indicated 
that the surface tension of the water droplet decreased drastically from 
72.7 dyne cm-"* - to less than 50 dyne cm"-*- at 20° C for droplets contaminated 
with grape drink and strawberry drink (Hung et al . , 1976). Therefore, 
the decrease in the surface tension of the water droplet exposed to a 
polluted atmosphere is taken into account. However, the order 
estimation of two terms in Equation (2-10) between a/r and e/r^ for the 
condition of a polluted atmosphere indicates that a/r is at least an 
order, or more, smaller than c/r . This means that the variation of the 
surface tension is not significant in the fog formation under polluted 
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Nevertheless, 


the expesut c vt the water di up lets 



atmospheric conditions, 
iu a polluted atmosphere makes a surface coveted with a muuolayei , h 
drastically decreases the effective difii.sicn cue ! i irieut , i leosUdi 
( 2 - 11 ), of the water droplet and inal.es the evaporation tat*.' slow 00,11. 


-f 

T 
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in. s'in>Y of tiii; NTMmoAi. limitation 


Tin 1 analysis at thr chemical compos it ion <0 atmospheric aerosols, 

the large uuantities ol St ), , NO., Nil! ami a strong acid :;ndi .if. H.hti, 

« a •* . a 

liavt* boon detected in highly indust rial i;s d ureas U.t ado re r et a!., 

1978). Strong evidence exists that tin* major producer oi No ha: • been 
tlu' automobile. There if. no doubt that mat ionurv si'ui.t lm! combust, ions 
also font r ilmt e a heavy traction ol tin* total NO, omit tod to tho utr. 10 - 
sphere. Fuel-bound nitrogen compounds an* ammonia, pyridine, and umv 
other amines (blassman, 19//). Similar to NO pn*duot ion by combust Inn, 

X 

sillier as an inhoront impurity in both ooal and petroleum, HO, if ! alst> 

A 

produfod through iuel oombust ion. 

There are .some very bu...c dill eivnee:. between the S problem and 
that post'd bv the lormutiou ol NO . Th * two possible :iuiin'i“. 01 N in 
any oombust ion profess are either atmospherlo or organ leal lv hound. 

S can be present, in elemental i'orm, organically bound, or as a f.peoief. 
in various inorganic compounds. Once it enter.*! tin* oombust ion process 
it is very reactive with oxidising species, and similar to fuel N, its 
conversion to the SO is fast compared to the other euergv releasing 
reactions (tflussman, 1977). The great quantity of atmospheric aerosols 
is thus formed from the pollutants of NO, and HO through the gas-phase 
photooxidation reaction. 

The average mass concentration of rosol particles from the National 
Air Sampling Network in the United States ranges from about 20 eg m ^ in 
clean air in remote areas, to values of b() to 2U0 sg m * in urban areas 
(Hidy and Brock, 1971). In the heavily polluted atmosphere of the highly 
industrialized areas, mass concentration with values as high as 2000 
pg m \ have been reported (He ini el d, 197')). As for particle number 
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* «= 'i , . 

concentrat ion, it varies i ram 10 particles * m in vvrv « jumii air t>> 

10* p art ivies urT * in polluted air. The ttdius of tenth umat*n uni 1<» 
with combust inn as a source ranges iron u.ul to 1 cm <!#» iaieid, i 9/b) 

The average particulates percentage distribution of tuns-, owniint '» 
tor SO is about 182, and for NO it is about (i (Hidy aim rOteb, Ib/lj. 

X 5* 

This im ans that mass concentration for 80 i;« roughly 1c;.. than 1.0 gg m 

in clean air, and ranges from about 20 to 90 eg i". ‘ in polluted ur. hie 

range oi mass concentration for NO. is roughly one* third el fin a;; un.* >0 

x 

SO . 

x 

In studying fog droplet form, t ion based on aerosol population, fear 
steps are considered. It can be shown as loll own: 

Step I. Muuodisperse, Single- Component Aerosol Model; 

Step II. Monodisperse, Multi-Component Aerosol Model; 

Step III. Nou~Monodisper.se, Single-Component Aerosol Model; 

Step IV. Non-Monod isper se , Multi -Component Aotvi ol Mulct. 

Step I. was adopted by Hung ot a l. (1979) and Hung and I. Jaw (19/9). Step 
II was considered by Hung and Li aw (1980). In the premtnt st.udv, Stet. ill 
will be described by Chapter IV, and Step IV will be uisot..,i>ea in then > 
V of this report. 

The numerical techniques of the application ot macroseop it Equations 
(11-1) to (2-9), with the inclusion of the evaluation of turbulent dJffut i 
parameters and other related physical parameters for the study of the dyn 
behavior of warm fog, are given in Hung and Vaughan (1.9/7). As lor the s. 
of numerical computation, the method of implicit integration oi the fin t 
difference equation was used to solve the governing paitial differentia! 
equations. In other words, the numerical technique employed is essential 
similar to that of Richtmyer (1957). The utilization of microscopic 
Equations (2-10)-(2-22) in solving the initial stage of t be nucleat ion 
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anti tlu* p.rowlh et t lit- water dreplets ,ihs.»c hit ed with a* r*«:,«*l parti* l* > i*- 
discussed In llunp, and llmklc ll'l/V); while tlu- prowl h and dot, n a? w.*f * t 
droplets in terms at seed i tip, «*i hvproscopie ehtmoals and a* t »*.•;*• I put i* St • 
ar»‘ dif.t ussed in llunp, i t al. UM/K). 

Result ly, llimp ft al. (1M?M) showed that tlu* aerosol port it lee v. i f Is 
tin* hip, her number density, larger sire nut lei, t lu* heavier mass as t he 
nuclei, and t hi* hip, bar rat ia at tin Van* t Ha-,' t ta* tar to tin mole. ular v.c ipht 
favor tin* formal ion at tin* lower visibil It v adveet ion tops with otronp.er 
Vertical enerpy transfer durinp thenueleed ian andcoudt n,,at ian t is;.* pat huh 

Tin* turbulent exchanp.e eo«*tt is lent.*. employed in this study are based 
at empirical t lux gradient relat ions, lit** detailed des.eript i.nts are listed 
in Appendix A. 

In this eatapntat ian, the method at implieit inieptation at the tinit. 
difference equal Ions has been ut ted to salve the poverninp partial di 1 1 - tent ia i 
equations. The initial values at the wind profile, temperature profile, and 
the relative humidity use tin* result.*; el field measurements taken a Ian**, t he 
California coast by Mack et al. (l‘t/1, 1 *J7 '♦ , 1 '»/•>)♦ The similar initial 
values of field measurements had also been used in tap, torra.it ion studies 
with detailed descriptions p.iven by llunp, et al. (1979). 

In the horizontal p.eomet rical X-eoordinate, the initial temperature 
at the p, round was 1 4 °C everywhere, and decreased to 12.3 0 C within inte and 
a half hours linearly at the location of the temperature sink 3 km * x • b km. 
The source terms of Equations (2-2) and (2-3) were assumed to be aero at 
time, t = 0; and I* - E at time, t U. Here, F = sF , where the 

CV CW CW CW-f 

subscript i denotes the species of hygroscopic aerosols servinp, as the 
condensation nucleus. F is computed from the time rate of clianpe of C! 

CW Wil 

defined in Equation (2-9) with the substitution of the p.rowth of a water 
droplet radius obtained from microphysical processes based on various 
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species of tiyp.ro scop ic aerosols given in Equation (2-lOj. Thus, the inicio- 
physics equation can be coupled with the macrophvsical dynamic equations 
through source and sink >.erms, F and F . Reasonably large time steps 
have been used to obtain a computationally stable solution. 

The visibility is computed by using the relation 

visibility » ™-|— (ra) (VI) 

where (■$ is related to the drop size distribution, N(r), and to Mle scattering 
efficiency, S^, which is important for droplets larger than 1 pm, according 
to 

00 

(5 = 4n S M N(r)rdr. (3-2) 

o 

S M is very close to the value unity if the drop size distribution is oi 
the same size. It takes a higher value, exceeding, perhaps 2, with a 
broad spectrum of the drop size distribution (Mason, 1971). To elucidate 
the fog optical properties, Low et a l. (1978) examined the effects of 
different spectral shapes and widths on gamma and lognormal distributions 
with the same liquid water content, and concluded that the drop size 
distribution greatly affects the result of visibility even though the liquid 
water content is the same. 

Fog is define in terms of visibility. Wnen the distribution of aero: '• i 
particles causes the visibility to go below the 1 km range the onset of fog 
has been determined (Silverman and Weinstein, 1974). In this study, the 
formation of fog implies that the visibility is below 1000m range. 

In this study, the radius of condensation nucleus, r, depends on 
the relative humidity S, mass of nucleus m^, and chemical composition 
in terms of Van't Hoff factor and molecular weight i/iq, , which may be 
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written (Mason, 1971) 


A. * t , a /r 

100 


(1 



(3 -3) 


whete a and c are expressed in Equations (2-15) and (2-16), respectively, 
In this study, computer program documentation, and Fortran listing, 
of the computer program are described and listed in Appendices 15, and 
C, respectively. 
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IV. ADVECTION FOG IN A POLLUTED ATMOSPHERE BASED ON 
MULTI -DISPERSE , SINGLE-COMPONENT MODEL 

In the present chapter we consider multi-disperse particle, size distri- 
bution as aerosol population in contrast to monodi.sper.se aerosol distri- 
bution considered by Hung and Liaw (1980). 

Three types of particle size distribution arc considered with keeping 
the total mass concentration as a fixed value. To simplify the calculation, 
and also to make the comparison of results easier, the combination of tluvc 
step functions for each distribution is chosen. These three types oi 
particle size distribution are expressed in figure 1. These distributions 
can be shown as follows: 

(1) Type A Distribution: 

N = 8100 particles/cm^, for r = 0.075 pm; 

N = 4000 particles/cm^, for r = 0.100 pm; 

N = 1700 particles/cm^, for r = 0.125 pm; 

(2) Type B Distribution: 

N = 3225 particles/cm^, for r = 0.075 pm; 

N = 4500 particles/cm^, for r = 0.100 pm; 

N = 2500 particles/cm^, for r = 0.125 pm; 

(3) Type C Distribution: 

3 

N = 2150 particles/cm , for r = 0.075 pm; 

N = 3000 particles/cm^, for r = 0.100 pm; 

3 1 

N = 3500 particles/cm , for r = 0.125 pm; 

where N stands for the number density, and r implies radius of nucleus. 

The total initial value mass concentration of the aerosol particles as 

condensation nuclei for Types A, B and C distribution are constant and 

3 

with a value of 54.9 pg/m . With keeping the mass concentration iixed, 
the total number of densities of Type A, B and C are 13,800, 10,225, 
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Figure 1. Three types of aerosol population distributions 
Each distribution is a combination of three ste 
functions in which the total mass concentration 
is conserved and equal to 54-9 gg/m- 5 . 
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NUMBFR DF-MSHY ( PART l Cl.ES/ cm j ) 



NUMBER DENSITY ( PARTICLES/ cm 3 ) 



NUMBER DENSITY ( PARTICLES/ cm 3 ) 



RADIUS OF CONDENSATION NUCLEUS DISTRIBUTION 
FIXED MASS CONCENTRATION = 54.9 *g/m 3 


and 8,650 partieles/em\ respectively. 

Since the pollutant, SO , dominates the combust ion-related air 

X 

pollution as indicated by Hidy and Brock (1971). The present study 
considers the contribution of SO as a major condensation nucleus tor 
fog formation, considering the effect of various non-monodispersu 
distribution aerosol population. 

In the present study, the initial values of the wind profile, tempera- 
ture profile, and the relative humidity use the results of field measure- 
ments taken along the California coast by Mack et al. (197.3, 1974 and 
1975) . The similar initial values of field measurements have also been 
used in fog formation studies, with detailed descriptions, given by Hung 
et al. (1979). 

In this study, the computation was carried out with one hour time 
coordinate after the fundamental equations were solved simultaneously 
in conjunction with the initial values adopted from field measurements 
provided by Mack et al. (1973, 1974 and 1975). The computation ol a 
one hour time coordinate was conducted for the purpose of making all 
the physical parameters involved in the model self-eonsistant. In the 
horizontal geometrical x-coordinate, the initial temperature at the ground 
was 14°C everywhere. After a one hour time coordinate computation, the 
ground temperature decreased one degree linearly in one hour at the 
location of the horizontal x-coordinate 3 km 5. x <, 6 km. The result of 
the time coordinate indicated in this paper, began at the time the ground 
temperature started to decrease linearly. 
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Comparison of the decree of fog formation in terms oi visibility is 
made for the variety of aerosol particle size distribution with total mass 
concentration of aerosols produced by energy-related combustion kept at a 
constant value. Ammonia sulfate is chosen in this paper as an example 
of the species of air pollutant SO which dominates the combust ion-related 
air pollution as indicated by Hidy and Brock (19 71) . 

Figure 2 illustrates relative humidity and visibility profiles at 
two horizontal locations, x = 3 and 0 km, and both aL vertical location, 
z = 1 m, of the formation of advection fog associated with type A distri- 
bution, shown in Figure 1, of the aerosol as condensation nuclei. Figure 
3 shows the similar relative humidity and visibility profiles of the 
formation of advection fog associated with type A distribution aerosols 
for different vertical location at z = 4 m, wiLh the same conditions as 
Figure 2. These two figures show that advection fog with visibility 
below 1000 m was formed at time t = 52 min. and relative humidity s = 

98.6% for altitude z =lm; and at time t = 86 min and relative humidity 
s = 97.8% for al'itude z = 4 m associated with type A aerosol particle 
size distribution. 

Similar computations are also accomplished for advection fog formation 
associated with aerosol particles, (NH^^SO^ as condensation nuclei with 
size distributions types B and C, as shown in Figure 1. Figures 4 and 5 
illustrate similar relative humidity and visibility profiles af two different 
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TYPE A DISTRIBUTION 
CNH 4 ) 2 S0 C| : N - 8100/cm 3 , R « 
N = 2/OO/cm 3 , R «» 
HASS CONCENTRATION = 50,9 

RELATIVE HUMIDITY — 


ALTITUDE 1 H 

,0/5 ; N n OOOU/cm 3 , R - 

,1/5 mm 

/i"/M 3 

VISIBILITY 


,100 




Figure 2. Visibility and relative humidity profiles of fog 

formation due to condensation nuclei, (NH,)„SO,, 

A 2 A 

with type A aerosol population distribution, as 
shown in Figure 1, at two horizontal locations, 
x = 3 and 6 km, and both at vertical location, 
z = 1m. 
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RELATIVE OIDITY 



TYPE A DISTRIBUTION ALTITUDE - R m 

(fIH/ 4 )oS0/. : N - 8100/cm 3 , R - ,075 ; N ~ ROOO/cm 3 , R = .lOO^M 

R = 1700/cm 3 , R - .125 mm 
MASS CQNCENTRAl ION » 5R.9 m^/m 3 
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Figure 3. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, (NlI^) 0 SO^, 

with type A aerosol population distribution, as 
shown in Figure 1, at two horizontal locations, 
x = 3 and 6 km, and both at vertical location, 
z = 4 m. 
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altitudes z = 1 and 4 m, respectively, with type I> aerosol size distri- 
bution, and the rest of the other conditions the same as that of Figure 2 
for the formation of advection fop, with ammonia sulfate as condensation 
nuclei. These two figures show that advection fop with vie. i hi J it/ below 
1000 m was formed at time t = 53 min and relative humidity s - 99. Uf' for 
altitude z - lm; and at time t » 92 min and relative humidity s 98.111 
for altitude z = 4 m associated with type B aerosol particle size distri- 
bution. Comparison between Figures 2 and 4, and Figures 3 and !j indicate 
that condensation nuclei with type A distribution initiates the formation 
of advection fog at least more than 1 min at altitude z = 1 m, and 6 min 
at altitude z = 4 m earlier than the condensation nuclei with type E 
distribution in this particular case. Furthermore, relative humidity 
at which fog starts to form is lower for aerosols as condensation nuclei 
with type A distribution than type B distribution. By re-examining Figur> 
1 to 5, it can be concluded that the effect of particle size distribution 
of aerosols as condensation nuclei is significant enough to affect the 
degree of fog formation, and also that the distribution with higher part, it 
number and small particle size is more favorable in dense fog formation 
than the distribution with lower particle number and larger particle size 
aerosols when the value of mass concentration is kept constant. 
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TYPE B DISTRIBUTION /[lllUlt = ■! h 

{NH.,) 2 SO„ ; A - 3?2Vwl , R ■ .071 m m ; N - 4W0/M 3 , R - .!;•!) m m 
H = 25G0/cm 3 , R - .125 mm 
MASS CONCENTRATION * 50,9 /.oV 
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Figure 5. Visibility ami relative humidity pro Files of fog 
formation due to condensation nuclei, (NIL ).,SO, , 

h «- t 

with type B aerosol population distribution, as 
shown in Figure 1, at two horizontal locations, 
x = 3 and 6 km, and both at vertical location, 
z = 4 m. 
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As to advent ion fog formation associated with tvpe t distr ifmt »»m 
of aerosols as condensation nuclei, Figures b and 1 .show similar relative 
humidity and visibility profiles hu two diiierent vertical Jurat i-ms, 

/. = 1 and 4 m, respectively, with the same conditions as Figure and i, 
except the diiierent particle sine distribution as indic.it id. These two 
figures a pa in show tltat advection fop with visifiilitv below lUUu was 
formed at t - min and relative humidity s = y«t. for altitude at 

z = 1 m; and at t ime t. - 44 min and relative humidity s c ho.'i lor 
altitude at z ~ 4 m. In comparison with Figures; 2, 4 and ft, the results 
show that (Condensation nuclei with type A distribution initiates the 
formation of advection fop, at least more than i.'j min; and type M distri- 
bution, at least more than i’,.3 min earlier than the condensation nuclei 
with type C distribution at altitude ;; ~ 1 m in this particular case, lor 
the altitude ■/. ■■ 4 m in comparison wiLh Figures 1, '> and /, the results 
show that condensation nuclei with type A distribution initiates the 
formation of advection fop, at least more than b min; and type 15 distri- 
bution, at least more than 2 min earlier than the condeusat ion nuclei 
wiLli type C distribution in this particular case. Furthermore, relative 
humidity at which fop starts to form is lower for aerosols as condensation 
nuclei witli type A distribution than type 15 distribution, which are in 
turn lower than type C distribution. 

Among the effects of air pollution on atmospheric properties, the 
most noticeable effect is perhaps the reduction in visibility, with and 
without the occurrence of condensation. The series of the present study 
is particularly interested in how aerosols, produced hv energy-related 
combustion, as condensation nuclei affect the formation oi advection 
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Figure (). 


Vi nihil iLy and relative Imui.litv proillos oi f.. - 
oimaUon due to condensation nuclei, C.i^ V-^u, , 

wxth typo c; aeronol population distr ilmt ,011" a! 

;: h l’ WI ! lu , I ' 1 « uri ‘ l. at UVO horizontal Iona 

* ~ , 3 aiul 6 kri » aiul both at vertical location 
7 . - i m. ’ 
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Figure 7. Visibility ami relative hnuidit y o 1 i 1 oi in 
formation due to comlensat ion mu* lei, (*,'H ,) .So, 

a .. a 

with type C. aerosol population distribution, as 
shown in Figure 1, at tv;o horizontal locations 
” i and ti hm, and both at veitieal location, 



fog. Hung and Liaw (1980) showud that Liu* greater the mass concentration, 
the particle concentration, and the size of Lhe radio:, of cond. .nation 
nuclei, the denser the fog is that is formed. However, by using the mono- 
dipserse distribution model, Hung and Liaw (1980) indicates that the particle 
concentration, rather than the size of the radius of nuclei provides 
greater contribution to the formation of fog, if the value of mass con- 
centration is kept constant. 

Observation shows that particle size distribution for combustion 
related aerosols are widely spread rather than monodisperse distribution. 

The results given in this study are based on multi-disperse distribution 
which is the combination of the three step functions for each distribution, 
and a single component aerosol model. The following conclusions resulted 
from the present study: 

(1) Aerosol distribution with higher particle concentration, rather 
than the size of aerosol nuclei, support a greater contribution 
to the formation of fog, if the value, of mass concentration is 
kept constant. 

(2) Relative humidity, at which advection fog with visibility below 

000 m is formed, is lower for aerosols of distribution with 
the higher particle concentration, which is in turn the smaller 
size of aerosol nuclei, than lower particle concentration, which 
is in turn the larger size aerosol nuclei, if the value of mass 
concentration is kept constant. 

(3) Aerosol population with high particle concentration-shifted 
distribution [see Figure 8(a)] provides more favorable con- 
ditions for the formation of dense fog than the aerosol 
population with low particle concentration-shifted distribution 
[see Figure 8(b) for description] . 
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V. ADVKCT10N FOG IN A POLLbTLl) ATMOSPHERE BASKI) UN A 
MULT I -DISPERSE, MULTI-COMPONENT MODEL 

In coni rant; Lo monodisperse aerosol distribution considered by Hun)’ 
and I,iaw (1980), the present chapter adopts multi -disperse particle size 
distribution as the aerosol population and the multi-component aerosol 
model to investigate the behavior of fog formation. 

For each different species, three types of particle size distributions 

are considered with keeping the total mass concentration of the aerosol 

particles as a constant value. To simplify the calculation, and also to 

make the comparison of the results easier, a combination of the three step 

functions for each distribution is chosen. Two components of pollutants, 

(NH^) 2 ^ 0 ^ and CaCNO^^* each component with three types of particle size 

distribution, are expressed in Figure 9. The total mass concentration 

for each component of pollutants is fixed. They are 54.9 and 13.3 pg/m"* 

for (NH^) 2 ^ 0 ^ and Ca(NO^) respectively. The range of mass concentration \ 

for CaCNO^),; is roughly one-third of the amount of (NU^^UCL. These 

distributions can be shown as follows: ' 

$ 

1 

(1) Type A Distribution \ 

i 

(NH^J^SO^ N = 8100 particles/cm^, for r = 0.075 pm j 

N = 4000 particles/cm^ , for r = 0.100 }im | 

3 | 

N = 1700 particles/cm , for r = 0.125 pm j 

CaCNO^),^ N = 1800 particles/cm^, for r 0.075 pm 

N = 1250 particles/cm"^, for r - 0.100 pm 

3 

N = 600 particles/cm , for r ~ 0.125 pm 


31 




(2) Type B Distribution 


(nh 4 ) 2 so 4 

N 

= 

3225 

particles/ cm*, 

for 

r -• 

0.075 

p:n 


N 

SC 

4500 

particles/cm*, 

for 

r - 

0.100 

pm 


N 

S3 

2500 

particles/cm*, 

for 

r = 

0.125 

]r.n 

Ca(N0 3 ) 2 

N 

SS 

1000 

particles/ cm"* , 

for 

r = 

0.075 

pm 


N 

= 

1200 

particles/ cm*. 

for 

r <= 

0.100 

pm 


N 

- 

800 

particles/ cm*, 

for 

r = 

0.125 

pm 

(3) Type C Distribution 






(nh 4 ) 2 so, i 

N 

E£ 

2150 

particles/cm* , 

for 

r = 

0.075 

pm 


N 

SS 

3000 

particles/cm , 

for 

r = 

0.100 

Pm 


N 

= 

3500 

particles/cm , 

for 

r = 

0.125 

pm 

Ca(S0 3 ) 2 

N 

SS 

550 

particles/ cm*, 

for 

r = 

0.075 

pm 


N 

= 

900 

particles/cm , 

for 

r = 

0.100 

pm 


N 

= 

1050 

particles/cm*, 

for 

r = 

0.125 

Pm 

where N stands for 

the number density, and r 

implies : 

radius 

of nucleus 


In Figure 9, the solid lines express the particle size distribution for 
(NH^^SO^; and the broken lines show the particle size distribution for 
CA(NO^) 2* With keeping mass concentration fixed, the total number 
densities of (NH^^SO^ for Types A, B, and C are 13,800, 10,225 and 8,630 
particles/cm , respectively; and the number densities of Ca(N0 3 ) 2 f° r 
Types A, B, and C are 3,650, 3,000, and 2,500 particles/cm , respectively. 

In the present study, the initial values of the wind profile, tempera- 
ture profile, and the relative humidity use the results of field measure- 
ments taken along the California coast by Mack et al. (1973, 1974 and 
1975) . Similar initial values of field measurements have also been used 
in fog formation studies, with detailed descriptions, given by Hung 
et al. (1979). 
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The numerical study of the application of macroscopic equation;;, 
a;; si town in Equations (2-1) - (2-y) , with the inclusion oi the evaluation 
of turbulent diffusion parameters and other related physical parametcis 
for the study of the dynamical behavior of warm fog, is given nv Hun;’ 
and Vauyhan (1977). As for the scheme of numerical computation, tin; 
method of implicit integration of the finite difference equation was; 
used to solve the yoverniny partial differential equations. In other 
words, the numerical technique employed is essentially similar to that 
of Richtmyer (1967). The utilization of microscopi; equations in 
solving the iniLial stage of the nueleation and the growth ol the water 
droplets associated with aerosol particles is discussed in Hung and 
Huckle (1977); while the growth and decay of water droplets in terms 
of seeding of hygroscopic chemicals and aerosol particles is discussed 
in living et al. (1978). Mathematical expressions of the computation of 
the visibility of fog from the drop size distribution of the waler drop 
spectrum, are given in Hung et al. (1979). 

For the purpose of making all the physical parameters involved in 
this model self-consistant , the computation of a one hour time coordinate 
was conducted after the fundamental equations were solved simultaneously 
in conjunction with the initial value adopted from field measurements 
provided by Mack et al. (1973, 1974 and 1975). In the horizontal 
geometrical x-coordinate, the initial temperature at the ground was 14°C 
everywhere. After a one hour time coordinate computation, the ground 
temperature decreased one degree linearly in one hour at the location 
of the horizontal x-coordinate 3 km <_ x £ 6 km. The results of the time- 
coordinate indicated in this paper, again, began at the time the ground 
temperature started to decrease linearly. 
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Comparison of the degree of fog formation in terms of visibility is 
made for a variety of aerosol particle size distribution for multi- 
component species with the total mass concentration lor each species of 
aerosols produced by energy-related combustion kept at a constant value, 
both ammonia sulfate and calcium nitrate are chosen as the example oi 
the species of air pollutants SO and NO which dominate the combust ion- 

X X 

related air pollution as indicated by Hidy and brock y 1 ‘ 1 7 ,L ) . 

Figure 10 illustrates relative humidity and visibility profiles at 
two horizontal locations, x = 3 and 6 km, and both at the vertical 
location z = 1 m, of the formation of advection fog associated with 
type A distribution, as shown in Figure 9 of the aerosol as condensation 
nuclei. Figure 11 shows a similar relative humidity and visibility pro- 
files of the formation of advection fog associated with type A distribution 

aerosols for different vertical locations at z = 4 m, with the same con- 

ditions as Figure 10. These two figures show that advection fog with 
visibility below 1000 m was formed at time t = 26 min and relative humidity 
S = 96.1% for altitude z = 1 m; and at time t = 61 min and relative 

humidity S = 95.6% for altitude z - 4 m associated with type A aerosol 

particle size distribution. 

Similar computations are also accomplished for advection for formation 
associated with two-component aerosols, (Nll^^SO^ and t'a(NO ^)., , as con- 
densation nuclei with size distributions of types 11 and C, as shown in 
Figure 9. Figures 12 and 13 illustrate similar relative humidity and 
visibility profiles at two different altitudes z = 1 and 4 m, respectively, 
with type B aerosol size distribution, and the rest of the other conditions 
the same as that of Figure 10 for the formation of advection fog with both 
ammonia sulfate and calcium nitrate as condensation nuclei. These two 
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TYPE A DISTRIBUTION 

(NH/j) oSOf, : N = 8100/ct^ , R - .075 an ; 

fl - «Q/cm* , R - ,1Q0 «m ; 

fl ° 1700 /ch 5 , R - .125 <-M ; 

MASS C'OHCENIPATION = W.^g/m 3 

RELATIVE HUMIDITY 


ALTITUDE - 1 m 

CA(tiOj) ? : N ° 1800/cm] , R •" /'A „m 
H = 1250/cm 3 > R " . 1UU «..m 
N " COO/cm 5 / P r . 17 j 
MASS CONCENTRATION « 18. > ,-u/f*. 3 

VISIBILITY ------- 



TIME ( min ) 



TIME ( min ) 


C*i 


Figure 10. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, (NH^f^SO^ 

and CaCNO^),, , with type A aerosol population 

distribution, as shown in Figure 9, at two 
horizontal locations, ;c = 3 and 6 km, and boLii 
at vertical location, z ~ lm. 
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TYPE A DISTRIBUTION 

(NtUoSOh : H - 8100/cm 3 , R = ,075 <im ; 

N - , R . ,100 «M ; 

N - 1700/cm 3 , R - .125 <_m ; 

HASS CONCENTRATION * M.9/ig/m 3 

RELATIVE HUMIDITY 


ALTITUDE - A m 

CA(N0 3 ) ? : • 1800/cm 3 , R - .0/5 //m 

' ' N " 1?50/cm, j R = .100 Hfi 
N - LOn/c-v' , R " 

MASS CONCENTRATION = 18. 3 W:v' 

VISIBILITY ------ 




TIME ( min ) 


Figure 11. Visibility and relative humidity profiles of fog 
formation due to condensation nuclei, (NH^SO^ 

and Ca(N0 3 ) 2 , with type A aerosol population 

distribution, as shown in Figure 9, at two 
horizontal locations, x = 3 and b km, and both 
at vertical location, z = 4m. 
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TYPE B DISTPI3UTI0N 

(;iHf|) 2 S0f i : li » 32? r -/t* 3 , R - .0/5 «m ; 

R .r -T;/C< , R - .lQfv.vi ; 

R - 2Y *9 , r < - .!?'• ; 

mss c-N* LN r i-\i : * r^.Q^v 


All I Till, - 1 m 


CACfJC 


V? 


N = iojo/cm ; , i: 

m - , k 


N - -/f •• - , • 

MAf'' I »|. f l(Tf AT 1 . ' . » < ; 

t\ t „»•> ». .t r in ’* - • * * » . 


Rf,LAT!Vr: HUMIDITY ■— 


VISIRIUIY 




Figure 12. Visibility and relative humidity profiles of log 
formation due to condensation nuclei, 

and Ca(NO^) 2 » with type B aerosol population 

distribution, as shown in Figure 9, at Lwo 
horizontal locations, x ~ 3 and 6 km, and both at 
vertical location, z = lm. 
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TYfT B DISTRIBUTION 

(Nil^SOj, : II - 3223 /cmJ , R - .073 mm ; 

N r: H’ffl/ct'Z , R * ,203 (IK ; 

N zKQ/cm^ , R - ,125 .k ; 

MAGS CONCENTRATION * 

RELATIVE HUMIDITY 


AIT I PILL Am 

CA(fiO*)o : M W!/'* 5 , R t- ,0/3 «, m 
N = l^'Yc-V , P * j' f iM 
N » 83 ,V4 3 , P. * .1:3 ■ m 
MASS C'JNCLNtPATIGN - 18,3 Wm 5 

VISIBILITY — ■ 



r 

1 

Figure 13. Visibility and relative humidity profiles of fog 

formation due to condensation nuclei, (Nil.) -SO. 

4 z 4 

and Ca(NO^) 2> with type B aerosol population j 

distribution, as shown in Figure 9, at two hori- 
zontal locations, x = 3 and 6 km, and both at 
vertical location, z ~ 4m. 

V , j 
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l ip, tires 


show that advoc l ion liv 


f Hi vi-ifu I it v helot; louft n wa- > i *» i s:s* *i 

at tine t =s Uuu win ami relative humility .*1 - b/,t» lor altitude / 1 m; 

and at tine t » #»7 inin and relative humidity S - *i »>. 2 lux altitude a 4 : 

associated with type 11 aerosol particle si/.e di. sir Shut ion. t.'oiupai i sous 

between Figure: 10 and Id tnd Figures 11 and 11 indicate that cmulen at i-*n 

nuclei with type A distribution imitates the 1 urination of .idvect ion lop 
at least more than 10 win at altitude z r -- l ra, and In win at altitude z 
Am earlier than the condensation nuclei with type 11 distrihutiou in 
this particular case. I’ur therm- re, relative humidity at which lop, starts 
to form is lower lor aerosols as condensation nuelei with type A distri" 
bulion than tvpe 11 distribution, by re**o:-:amiiiing Figures 1-3, it can 
be concluded that the effect of the particle sine distribution ol t lie 
aerosols as condensation nuelei is significant enough to uliecL the dopjer 
of fop formation and also that the distribution with a sip, her particle 
number and a smaller particle size is more favorable in dense fop lormat ion 
than distribution with a Lower particle number and larger particle sice 
aerosols when the value of mass concentration is kept constant, ibis 
result further coni irms the computation obtained lor the various aetosol 
size distributions wiLh the s ingle-component model. Comparisons between 
s i.np,l e - component and multi -component models show that the mult i -component 
model definitely enhances the dense lop, formation more than the siup.le- 


component model. 

As to advuction fop, formation associated with type C distribution 
of aerosols as condensation nuclei, Figures 14 and 13 show similar relative 
humidity and visibility profiles for two different vertical locations, 
z-1 and 4 m, respectively, with the same conditions as Figures 10 and 11 
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Figure 14. Visibility and relative humidity proiiles oi i 
formation due to condensation nuclei, (Nii^).,Ho 

and Ca(N 0 ,j> 2 » with type C aerosol population 

distribution, as shown in Figure 9, at, t wo 
horizontal locations, x " '3 and 6 kin, and both 
at vertical location, ?. ~ 1m. 
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with the exception of the different particle size distribution an indi- 
cated. These two lip, ores again show that advection foe, wiLh visibility 
below 1000 m was formed at t = 38 min and relative humidity S = 97.22 for 
altitude at z * Ira; and at time t » 71 min and relative humidity 
S ® 96.42 for altitude at z = 4 in. In comparison with Figures 10, 12 
and 14, the results show that condensation nuclei with type A distribution 
initiates the formation of advection fog at least, more than 12 min; and 
type 15 distribution, at least more than 2 min earlier than the condo 
nation nuclei with type C distribution at altitude z = 1 m, in this 
particular case. For the altitude z - 4 ra, in comparison with Figures 
11, 13 and 15, the results show that condensation nuclei with type A distri- 
bution initiates the formation of advection fog at least more than 20 
min, and type B distribution, at least more than 4 min earlier than the 
condensation nuclei with type C distribution in this particular case. 
Furthermore, the relative humidity at which fog starts to form is lower 
for aerosols as condensation nuclei with type A distribution, than type 
B distribution, which in turn, is lower than type C distribution. 

The results given in this chapter are obtained from multi-dxsperso 
distribution, which is the combination of three step functions for each 
distribution, and a multi-component aerosol model. The following con- 
clusions resulted: 

(1) Aerosol distribution with a higher number density, rather than 
the greater size of aerosol nuclei, produce a more favorable 
contribution to the formation of lower visibility of fog, if t de- 
value of mass concentration of aerosols is kept constant. 

(2) Relative humidity, at which advection fog with visibility below 
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1000 m is formed, is lower for aerosols of distribution with 
a higher particle concentration if the value of mass; concent rut ion 
of aerosols is kept constant. 

(3) In addition to the contribution made by SO which dominates the 
air pollution produced by energy-related combustion, No , which 
is roughly one-third of the. amount of SO , provides a considerable 

X 

contribution to fog formation, and is important enough to lie 
included in the fog formation prediction. 
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VI, LIQUID WATER CONTENT RESULTING FROM CONDENSATION EFFECTS 


The mathematical simulation of the formation of advection fog based 
on the multi-disperse distt Lion of the aerosol population shows that 
the aerosol distribution with a higher number density, rather than greater 
size of aerosol nuclei, produce more favorable conditions for the formation 
of denser fog, if the value of mass concentration of the aerosol particles 
is kept constant. It is of considerable interest to study the liquid 
water content resulting from the effects of the condensation processes in 
which the aerosol particles with a variety of distribution functions are 
taken into consideration. The results of liquid water content are very 
important in the study of the acid rain problem. 

A growing body of evidence suggests that acid rain may have sub- 
stantial adverse effects on the environment. Such effects include acidi- 
fication of lakes, rivers, and groundwaters, with resultant damage to 
fish and other components of the aquatic ecosystem; acidification and 
demineralization of soils; reduction of forest productivity; damage to 
crops; and deterioration of man-made materials. These effects may be 
cumulative or may result from peak acidity episodes. The sLudy of liquid 

water content due to condensation processes associated with SO and NO 

x x 

as condensation nuclei of multi-disperse distribution may be beneficial 
to the study of potentially impacted ecosystems due to acid rain. 

With the same initial and boundary conditions we considered in the 
previous section, the cime dependent liquid water contents were calculated 
with types A, B and C particle size distributions, as shown in Figure 9 
for (NH^J^SO^ and CaCNO^^ as two major components of the aerosol particles 
as condensation nuclei. Figure 16 illustrates the time change of the 
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Figure 16. Time rate of change of the amount of liquid water content 

associated with condensation nuclei, (NH,) SO. and Ca(NO ) 

■4 1 4 3 

with type A aerosol population distribution, as shown in 
Figure 9, at two horizontal locations, x = 3 and 6 km, and 
both at vertical location, z = 1 m. 
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LIQUID WATER CONTENT ( g/m 3 ) 



TYPE A DISTRIBUTION ALTITUDE -1m 

(NHij^SOjj : N = 8100 /cm 3 , R = ,075 >/jm; CACNO?^ : N = 1800/cm 3 , R = . 075 u m 

N = AOOO/cm 3 , R = .IGO/zm; N = 1250/cm 3 , R = .100//M 

N = 1700/cm 3 , R = . 125/1 m; N = 600/cm 3 , R = ,125/im 

NASS CONCENTRATION = 5A.9 //g/m 3 MASS CONCENTRATION = 18.3/<g/m 3 


increase of the liquid water contents at two horizontal locations, x : j 
and 6 km, and both at the vertical location z - 1 m, of the formation ut 
atmospheric acid deposition of wet removal processes associated with type 
A distribution of aerosol as condensation nuclei. Figure 17 shows a 
similar time change of the increase of the liquid water content of tin- 
formation of acidic deposition of wet removal processes associated with 
type A distribution aerosols for different vertical locations at z - 4 m, 
with the same condition as Figure 16. These two figures show that acidic 
deposition of the liquid water content becomes 102 yg/m^ at time t = 40 

3 

min for altitude z = 1 m; and 130 yg/m at time t « 75 min for altitude 
z = 4m associated with type A aerosol particle size distribution. 

As to the type B aerosol size distribution, Figures lb and 19 
Illustrate a similar time change of increase of acidic deposition of 
liquid water contents at two different altitude z ~1 and 4 m, respect ively, 
with the remaining conditions the same as that of Figure 16. A comparison 
between Figures 16 and 18, and Figures 17 and 19 indicate that the formation 
of atmospheric acid deposition of wet removal makes no difference lor 
condensation nuclei with types A and B distributions. 

As to the time change of increase of liquid water content associated 
with type C distribution, Figures 20 and 21 show similar profiles at two 
different altitudes, z =1 and 4 m, respectively. Comparison with Figures 
16, 18 and 20, and Figures 17, 19 and 21 show that the atmospheric acid 
deposition for wet removal processes make no difference for the aerosol 
population with any distributions so long as the mass concentration of 
the aerosol particles is the same value. 


47 


£ T. T, 

in cd ia 

o cn 

O rH rH 

II II II 
c< oi OS 


rA n-\ r/~\ 
K K X 
u u u 
\ ^ 
n a o 
o in o 

00 CM UD 


K 

u* 

ii 

LU 

CD 

ro 

h— 

— j 
< c 


II I! II 


CN 

(A 

O 


<C 

CD 


5 : 5 : sn 

^ =C V 

LA CD LA 
l> CD CN1 




o 

rH 

> — i 

CD 

- 



f— H 

h~ 

li 

II 

II 

CD 

PC 

ns 

CtS 

DC 

s 


■* 


1 — 

ND N"\ NA 

CO 

JT 


s: 

i — ■ 

U 

«D 

u 

CD 

N 

N 

N 


CD 

CD 

CD 

«r 

CD 

CD 

O 


» — 1 

O 


LU 

OO 

■=r 

rH 

a. 

>- 

II 

II 

II 






O 

5 


dT 

c n 

CN 

•=r 


ZT 

LA 


II 

O 


<c 

as 


o 

CD 

co 

co 

<c 


5 : 



2 

5 


UJ 

>t 


( £ W /9 J 


1N31N03 H3iVM aiflOH 



* 




CJ 






as 

4-J 

CA 


s 

q 

o 

r*< 


o 


•*H 


4-J 

v-/ 



q 

q 

i—* 


c 

u 

*> 


a 


J 



T3 


'w 

q 

q 

X 


a 



ad 

U 


X 

q 

q 


q 

q 

£ 






r« 



C4 

q 


*H 

/— s 

c 

il 

q 




cr 


4—1 


•H 




rH 

w 

rC 

r 



•H 

X 

cm 

* 

q 

q 

O 

♦H 

4-i 

a 


<y 

X 

*r-< 

4J 

rH 

•H 

4-J 

q 

u 

AD 

q 

q 

q 


o 

o 

q 

q 

o 

£ 


c 

1 — 1 

q 

q 

*r-l 



<u 

-c 

w 


o 

•H 

AJ 


CO 

c 

<y 

T3 

o q 
eo o 


o 


c 

q 

*c 

o 

U-H 

o 


O Z) 
4-J 4-» 

q q 
q *h 
u 

a o 

e x 

»H X 

H q 


cu 

q 

q 

oc 

•H 

fr* 


4D 

q 


c 

X 

c 

q 

O 


C- O'. 


q 

o 

»H 

4-1 

q 

o 

o 


q 

u 


q 

o 

> 


4-i o q 
q 


•-« — o 
2 ^ ^ 


48 


«E ST E 

LO C3 LA E 

o cnj <: 

O i — I i — | o 

’ ‘ ^ 

II II II 

az cc oo 


N“\ fv\ m 
2: s: 2: 
<-> u u 

N \ S. 
OOO 
o a o 
o CN 00 


y cm 

r 3 m 

H- CD 


— J -cc 

«=C CJ 



2; 5: 2: 

, woui s: 

^ ^ g ts a 

O O r-H r-H 3 


( £ W/9 > 1N31N0D H3JLVM QinCH 


h— 

=3 

II 

II 

II 

CD 

OQ 


cr 

cd 

QO 

cel. 

ac 

Ln 

1 — 

■% 

_ 

_ 

II 

co 

n rn ir\ 

►— < 


21 

2*1 

sr 

Q 

0 

<J 

U 

0 

.OQ 

s 

M 

\ 


lh 

O 

O 

H- 


CM 

CD 

CD 

<C 

UJ 

CM 

LT\ 

Ln 

QC 

Q.. 


■sr 

CM 

h-~ 


II 

II 

II 

UJ 





O 


z: 

sr 

0T 

or 





0 


•• 



0 


cr 



CO 


CD 



CO 


00 



<c 


CM 



2: 



CO 

o 'o 


aj m 
3 O 
3 w 

AJ 

C 3 

O CJ 

u 

T) 

>•< s 

3 3 

4J 

3 <j- 

S O 

C /3 

T3 Cvl 
•H * 

3 -J- • 


h> r ~ < 

D* BJ 

AJ 

K 


1 

•H S3 

3 



O 

rH 




r-i 


•H 

« 


Ln 

4-4 n 


d 


O *H 

AJ 

0 

» 


3 

W 

•H 

S 


Al M 

*H 

AJ 


c 0 

T 3 

d 

rH 


3 3 


u 



2 c 

d 

0 

11 


6 

O 

pH 


3 3 

•H 


N 


O 

AJ 

pH 



3 H 

d 

d 

r* 

. 

Xi AJ 

pH 

AJ 

d 


AJ 3 

d 

d 

0 


w wn -w 

ft'H # 
s-i o 
H O O 
O *3 r—4 
« 

O O r-4 

^ S 3 
3 aj o 
3 -H 

AJ A) 
(3 3)-i 
3 

3 * > 

O. O': 

>) 4J 
AJ 3 3 
l-i 

-3 3 JS 
A-> tO AJ 
•r 4 *H O 

& fa JO 


£ w / 9 ) 1N31N03 H31WI amon 


49 


TYPE B DISTRIBUTION ALTITUDE = H a 

(rJH i| ) 2 S0 £l : N = 3225/cm 5 R = ,G75jah; CA(N0 3 ) 2 : N = 1003/cm 5 , R = ,375/»m 

N = 45G0/cm 5 R = .100/iMj N = 1203 /cm 5 , R = .1S0«« 

N = 2500/cm , R = .125/imj N = 8 jO/cm 5 , R = ,125/<m 

MASS CONCENTRATION = 5A.9/<3/m 5 MASS CONCENTRATION = 18.3/*s/m 3 
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with Tvpe C aero- „>1 population distribution, as shown in 
Figure Q , at two horizontal locations, x = 3 and 6 Isa, and 
both at vertical location z = 1 m. 
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VII. CONCLUSIONS AND DISCUSSION 


The analysis of the chemical composition of atmospheric aerosols, 
the large quantities of SO^, NO^, NH^ and a strong acid such as H^CO^ 
have been detected in highly industrialized areas (Leaderer et al. , 
1978). Strong evidence exists that the major producer of NO has been 

X 

the automobile. There is no doubt that stationary source fuel combus- 
tions also contribute a heavy fraction of the total NO^ emitted to the • 
atmosphere. Fuel-bound nitrogen compounds are ammonia, pyridine, and 
many other amines ( Glassman, 1977). Similar to NO production by 

X 

combustion, sulfur as an inherent impurity in both coal and petroleum, 
SO^ is also produced through fuel combustion. 

There are some very basic difference between the sulfur problem 
and that posed by the formation of the nitrogen oxides. The two 
possible sources of nitrogen in any combustion process are either 
atmospheric or organically bound. Sulfur can be present in elemental 
form, organically bound, or as a species in various inorganic compounds. 
Once it enters the combustion process it is very reactive with oxidiz- 
ing species, and similar to fuel nitrogen, its conversion to the 
sulfurous oxides is fast compared to the other energy releasing reaction 
reactions (Glassman, 1977). The great quantity of atmospheric aerosols 
is thus formed from the pollutants of NO^ and SO^ through the gas-phase 
photooxidation reaction. 
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Among the effects of nir pollution on atmospheric properties, 
the most noticeable effect is perhaps the reduction in visibility, with 
and without the occurrence of condensation. The series of present 
study, which we have accomplished here, have been particularly 
interested in how aerosols, produced by energy-related combustion, as 
condensation nuclei affect the formation of adveetion fog. These 
aerosols provide the major source of condensation nuclei for the forma- 
tion of fog in the areas of busy highways, airports, seaports, etc. 

How these pollutants affect our living environment, in particular the 
formation of adveetion fog, is also our major concern in the present 
study. 

The conclusions obtained from the present series of studies can be 
summarized as follows: 

(1) Hygroscopic chemicals with higher ratio of the Van't Hoff 
factor to molecular weight produce more favorable conditions 
for denser adveetion fog than chemicals with lower values of 
the ratio 

(2) A higher number density of aerosol nuclei produces adveetion 
fog with lower visibility during the nurl cation and condensa- 
tion processes than a lower nuclei number density. 

(3) Heavier mass nuclei aerosol particles favor the formation of 
adveetion fog with a lower visibility than do lighter mass 
nuclei during the nucleation processes. 

(4) Condensation nuclei with larger radii cause a denser adveetion 
fog than condensation nuclei with smaller radii. 

(5) For a condensation nucleus to grow into a droplet, the air 
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must have attained a certain decree of supersaturation. This 
condition is not necessarily true for a condensation nucleus 
associated with a polluted atmosphere which can be easily 
grown into a droplet and produce a dense fog without having 
the air attaining supersaturation. 

(6) In the analysis of air pollutants, the range of mass concen- 
tration for NO is roughly one-third of the amount of SO . 

The major contribution of combustion-related pollutants as 
condensation nuclei in the formation of fog comes from 
aerosols with species of SO . However, NO as condensation 

XX * 

nuclei also gives a considerable great contribution to the 
formation of fog. 

(7) Aerosol distribution with higher particle concentration, 
rather than the size of aerosol nuclei, support a greater 
contribution to the formation of fog, if the mass concentra- 
tion is kept constant. 

(8) Relative humidity, at which advection fog with visibility 
below 1000 m is formed, is lower for aerosols of distribution 
with the higher particle concentration, which is .in turn the 
smaller size aerosol nuclei, if the value of mass concentration 
is kept constant. 

(9) Aerosol population with high particle concentration-shifted 
distribution provides more, favorable conditions for the forma- 
tion of dense fog than the aerosol population with low particle 
particle concentration-shifted distribution. 

The results of the present study are in reasonably good agreement 




with field observat ions from the California eoaot made by Ma< \ ot al . 
(1973, 1974 and 1975), and others, although a detailed ro.-.par 1' on if. 
not possible .it this time because of the Unit number of t tint* dependent 
measurements of dynamic renditions and aerosol data in the horizontal 
and vertical coordinates. 

The present study can be developed as a numerical predict ion oi 
the formation of advert ion fog with the 1 vale . of following parameter!; 
as initial input: (1) wind profile; (2) temperature profile; 

(3) humidity profile; (4) mass eoneent rat ion of aerosol part icle;.; 

(5) part icle sice distribution of aerosols; and (6) chemical composi- 


tion of aerosols. 
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APPENDICES 


APPENDIX A 


EVALUATION OF KEY PARAMETERS IN THE 
MATHEMATICAL MODEL 


The turbulent exchange coefficients proposed in this model are 
based on empirical flux gradient relations which have been measured in 
studies of turbulent transfer in the lowest 50 m or so of the atmo- 
sphere (Businger et al., 1971; Panofsky, 1974; Pasquill, 1974; Huang, 
1975). 

(i) Parameterized Fluxes of Momentum, Heat and Moisture 
The parameterized fluxes through the surface boundary layer serve 
as the boundary conditions for the turbulent exchange processes in the 
transition layer. The "bulk aerodynamic" method (Roll, 1965) was used 
to parameterize the surface boundary layer fluxes of momentum, heat 
and moisture. 
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Here u = (x/p) is the friction velocity; H, >:• .■? vertical heat 
flux through the surface boundary layer; E, the flux of vapor through 
the surface boundary layer; and T, the surface stress tensor. Within 
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the constant flux layer, it can be shown that the non-dimensional wind 
shear, temperature gradient and vapor flux gradient can be expressed 
in the form 
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scaling water vapor mixing ratio. 


By definition, then, within the surface boundary layer, the 
exchange coefficients for momentum, heat and diffusivity can be shown 
to be: 
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In the advection fog model, these turbulent exchange coefficients 

are multiplied by the factor exp [-8fz/u*] , obtained by Shir (1973) 

in a numerical study of a neutral boundary layer, to provide a 

reasonable behavior of the turbulent exchange coefficients above the 

constant flux layer in the model. 

In this model, the relations of non-dimensional wind shear 
z z 

4> (~f -); temperature gradient, <f>, (~r-); and water vapor flux gradient, 
in l n Li 

S s 

); are used. These parameters are based upon atmospheric 

s 

measurement within the surface boundary layer under relatively steady- 
state and horizontally homogeneous conditions ( Pancfsky, 1974) . 

Under stable conditions L g > 0, the relations (Webb, 1970) 
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are used. Under unstable conditions or I, <0, the relations (Dyer, 
1974) 
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are used. Under neutral conditions or H = 0. both the stable and 
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unstable relations reduce to (Dyer, 1974; Huang, 1979) 


*h - *i ■ ■ 1 <*-»> 

The turbulent exchange coefficients K^, and can then be 
computed for every grid point in the model from Equations (A-l) to 
(A-12), using the friction velocity u* and the scaling length L 
determined for the column. 

(ii) Infrared Radiation Flux Divergence 

The model for the infrared radiative flux divergence produced by 
fog adopted the model proposed by Goody (1964) 
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where cr is the Stefan-Boltzmann constant; K , spectrally-averaged 
s w 

2 

mass absorption coefficient of fog for infrared radiation (cm /g); 
and z rj ,» the height of the fog top. The absorption coefficient K is 
the model proposed by Zdurkowski and Nielsen (1969) 
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for the absorption per unit path length in the fog. 

(iii Terminal Velocity of Fog Drops 

The sedimentation of the fog drops can be simulated through the 
mean terminal velocity U fc which follows Stokes relationship (Fletcher, 
1966) 
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where r is the radius of the fog drops. 


APPENDIX 15 


PROGRAM DOCUMENTATION 

15.1. Introduction 

The program documentation provided in this Appendix is separated 
into several parts. (1) Within the Introduction the essential features 
of the program are described together with a discussion of some model 
capabilities, (2) A brief discussion of the control which can be 
exercised over the program in order to study particular features of 
advection fog. (3) A description of input cards. (4) A listing and 
definition of variables and controls used in the program. (5) Three 
pages of sample output. This presentation is designated to allow 
others to make the program accomplish their desired aims with a minimum 
expenditure of learning time. 

The program described here provides a tool for studying both the 
formation and dissipation of advection fogs in the vertical (X-Z) plane. 
In the discussion which follows* quantities which occur in the computer 
program are capitalized; those that are input quantities are underlined. 

The model can describe the evolution of the fields of U, V, Y, PT, 

R and W from a variety of initial conditions, subject to the imposed 
upper and lower boundary conditions. The initial conditions are put 
into the program by sets of input cards which are based on the field 
observations. The lower boundary conditions represent the earth's 
surface, either land or ocean, at which PT is constrained to be equal 


65 


to the surface, temperature. The uluee temperature can tail- on many 
desired configurations along the X-axis. At the upward boundary, the 
dependent variables are maintained equal to their Initial values. Al 
the downwind boundary, however, the dependent variables are continually 
adjusted to equal their computed value in the adjacent upwind grid 
column. 

Wihtin each time step (1JT) , the U, V, Y, PT, W and R fields are 
dumped by the simulated processes, in accordance with the boundary 
conditions, by an implicit integration scheme. The temperature (i ) 
field is then diagnosed from the potential temperature (PT) field. 

The occurrence of any condensation or evaporation (FG) is accounted 
for by the local conditions and is calculated by the microphysics 
equation subroutine RKF. From these final fields, the diagnostic 
variables, turbulent exchange coefficient (KA and KM), total water 
below a height (Z) in a column (TNT), specific heat of saturated air 
(CPT) , radius of droplet (RAD) and visibility (VIS1TY) are computed 
and stored for use in taking the next time step. 

Output is controlled by output time (OT) and ending Lime (liT) , 
which occurs when the total time elapsed (TIME) at the end of a set 
of time steps equals the print time (PRT) . Printout of the fields 
of wind velocity (U, V and Y), vertical eddy heat flux (HO), T, R, 
liquid water content, turbulent exchanges coefficient, radiative flux, 
radiative cooling rate, water vapor flux, radius of droplet and visibi- 
lity, in that order, is provided. Integration then continues, with 
output occurring at PRT equal to integer multiples of output time (OT) , 
until program termination at the ending time (KT). 
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15.2. Control of Program 

Flexibility of the program to investigate various aspects of 

•» 

advection fog formation and dissipation rests with the way the con- 
figuration of the program can be changed through specification of 
input data and control indices. These quantities fall into the follow- 
ing categories: 

a. time parameters 

b. grid spacing and configuration in the vertical and horizontal 

c. input parameters controlling which physical processes operate 
and their magnitude, and the initial values of the dependent 
variables 

d. boundary conditions on temperature and water vapor at the 
earth's surface boundary 

e. options for printout 


a . Time 

The basic time varibale is the length of the time step DT, which 
can be changed to insure computational stability, and control running 
time on the computer. Additional time variables are OT, for multiples 
of which the program prints out and ET, which is the total meteorolo- 
gical time the program is to run. 

b. Grid Spacing 

The vertical grid is composed of (KE) levels (up to a maxium of 
50) with the space between adjacent levels expanding based on the 
height of the first level ( ZAL ) above the earth's surface and an expan- 
sion factor (ZAK) . 

The horizontal grid is composed of ( IE ) columns (up to a maximum 
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of 40). JN£ can be specified as any number greater than two a»d < 40, 
except that printout Is conveniently grouped in units of ten columns 
per page. The horizontal grid is basically specified by ( UK LX ) which 
is the uniform spacing between columns in that portion of the grid 
located between ( IL) and ( IJt ) , the (I)-values of the left and right 
boundaries of the region of uniform grid spacing. Outside these limits, 
the grid spacing expands to the upwind and downwind grid limits with 
the size of spacing controlled by the expansion factor (XA1) and 01. LX. 

In the output, negative X's Indicate columns located upwind of the uni- 
formly spaced grid and positive values label the uniform portion and 
columns in the downward expanding portion, 
c. Initial Values 

If I SKI) ■ 1 and 1KAU * l, sedimentation and radiational cooling, 
respectively, are included in the model; if 1S1-.U and IRAQ i 1, then 
the processes are excluded. 

PT, W and R can be Initialized to be uniform everywhere having 
values ( PTI ) , (ZKRO) , and (RI_) , when the corresponding control indices 
( IPT , IW and IRR ) ■ o. If IPT ■ -1, then the temperature profile is 
isothermal with a value of (TP) and the correspondxng initial PT is 
computed from the hydrostatic relationship. If JTT, IKK , J_W and J_U 
all * 1, the corresponding varibales are uniform in X, and T, R, W, 

U and V are read from a card listing for K * I, .... KK; the prognostic 
variables PT then be computed. If JJJ =* 0, then U and V is an adiabatic 
profile at all columns, using friction velocity (UP) and roughness 
length (Z0). 
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of 40). IK can In* specified as any number greater than twc .cd < 411, 
except; that printout is conveniently grouped in units of ten eolumns 
per page . The horizontal grid is basically specified by (DKhX) which 
is the uniform spacing between columns in that portion of the grid 
located between (II.) and (IJR), the (I) -values of tlie left and right 
boundaries of t'he region of uniform grid spacing,. tJutsi.de these limits, 
the grid spacing expands to the upwind and downwind grid limits with 
the size of spacing controlled by the expansion faetor (XAI) and DI.I.X. 

In the output, negative X's indicate columns located upwind of the uni- 
formly spaced grid and positive values label the uniform portion and 
columns in the downward expanding portion, 
c. Initial Values 

If I SKI) - 1 and JRAD - 1, sedimentation and radiational cooling, 
respectively, are included in the model; if I_KKb and I HAD i- 1, then 
the processus are excluded. 

PT, W and R can be initialized to be uniform everywhere having 
values (PTI) , (XKRO) , and (RI) , when the corresponding control, indices 
( IPT , I_W and IRR ) ~ o. If IPT = -1, then the temperature profile is 
isothermal with a value of (TP) and the corresponding initial PT is 
computed from the hydrostatic relationship. If IPT, IRR , IW and 1JJ 
all = 1, the corresponding varibalos are uniform in X, and T, R, W, 

U and V are read from a card listing for K = 1, ..., KK; the prognostic 
variables PT then he computed. If _IU = 0, then U and V is an adiabatic 
profile at ail columns, using friction velocity (UP) and roughness 
length (Z0) . 
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il. Boundary Values 

The boundary rendition of potential temperature on the surface is 
controlled by the index (IimiM). If IOTKM - u, then I)TFM(i), the 
surface temperature change being investigated (negative tor adveetion 
fop, formation, positive for fop, dissipation) equal to zero, if 
1DTEM * 1, then DTKM(I) is uniformly equal to IffhMI between IT KMI. and 
TTKMR, tlie left and rigrt I- boundaries of this warmer or colder rep, Ion. 
If Jilf/iM “ -1, then DTKM(I) is read from a card list tor all i's except 
1*1 and IK. In all these cases, I)TKM(I) is introduced gradually to 
avoid pulsinp the integration. The length of time over which Dl'KM(I) 
is introduced is equal to (TIM) . 

If (I RS PC ) = 0, then a boundary condition of no vertical flex is 
imposed on R at the surface. If ( IRS PC ) * 1, then the boundary 
condition of R is set as saturation mixinp ratio at surface temperature 
and pressure equal to 1000 mb. If ( IRS PC' ) = 2, thou the boundary 
condition of R at surface is kept the same as the initial. 

Conditions at the upstream boundary remain unchanged throuphout 
the integration, while conditions at the downstream boundary are set 
equal to those at the adjacent column at the end of each time step, 
e. Output 

Output from the program is controlled by the main program and 
performed by a print subroutine (PRNT) . For variables in storage, a 
call of PRNT, with proper arguments, accomplishes the output. For 
variables not in storage, the values are computed in the main program 
and stored in a dummy storage location called HC, and output occurs by 
calling PRNT with HC as an argument. 
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flu* prim subroutine prints -» tl.it .t col win-. per pa-t- lor each 
output quantity, along wits documentary vertical ami horizontal * i 1 . . = 
tame information. I.ucii column rover;; tin* entire vertical extent oi 
the grid. bogie in the program provides for print ing only li, column. oi 
data when if is less than at) and an integer multiple of 111, In addit i<»u, 
the first page identifies the output variable, output time, and other 
selected informat ion which identifies the simulation. 

f. Units 

Internally, the program operates in cgs units, except that Ul !.X 
is in meters and when heat is an explicit mil it is expressed in 
calories. 

g . Computational Requirements 

The program has been written in Fortran and run on a Univue I lift 
digital computer. Approximately 20 K words of core storage were 
required for execution of the program. Cl’U time requiremeuis lor exe- 
cution of the program depend upon the ratio KT/DT, the total number 
of grid points (IK) x IKF.) , and the amount of output specified. In 
addition, more, computations per time step are required in fog than in 
non -fop, regions. 

The primary limitation on the maximum time step (I)T) that should 
be employed is the computational stability criterion, 

U * 1)T 

IJI-LX 1 ’ ( 1 ) 

for the integration of the horizontal advection terms. It is advisable, 
however, to restrict l)Tv.60 seconds, even when condition (1) can In- 
satisfied by a larger I)T, in order to avoid serious truncation errors 
and/or computational instabilities in the vertical integrations. 
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Ii . ,'i . Sample i_Sct of Input Cards 

a. lime Variable Constants, Temperarure and Wind Input 


0*7=-t *2 0.90K + 3 7 * 2 G£. + 3 7.2C£*3 ?.*«£-! «7i£ + l 2.07E-*' t.BOt-7- 


DT 

1091514 

1 1 < 1 1 i 

1 u i n i 


or 

99 199 

i 1 <i 1 t i , 
>111111 


KT 
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■I » I ‘ I'l 
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Ml 939 

, : ;i II I ,1 
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R r - 

9 6 9 (394 
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I I I II I 


KW 
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i ii ii ii ii i ii 

III II 


uuinnuin 


rr 

9 9 9 9 9 9 9 

ii v ii ii 11 ii 

II I II I I 

in:: 


A) 

99 


• 0 Of *1 S* 0 Cf < < 


8 0 . ft 9 0 1 9 

* 

1 1 1 1 1 ip i : 1 1 1 1 1 


Vi 

99 10 

nil' i i «i 

1 I I I I I I 


nnn : nnnnnnn 


1 1111113 1 1 1 3 33 3311 33 3333 11 1131131 1355333 ’313153 33333)5 113311) 33)1) 
ltllllllltOIIIUIi>l|i|UIIIII|IHOIIMUU(IMOIIIUIUUHUOUIIII|l|lll 
J J S 5 5 5 S 5 5 S 5 5 5 5 5 1 3 3 5 5 5 5 3 5 5 5 5 5 5 i 5 5 5 5 5 5 5 5 5 5 5 5 5 5 ; 5 , 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 
M ( 5 1 5 1 M 5 6 5 5 5 1 1 i ( C 6 I 6 ( t I I M S 5 5 5 f t ( 5 ( 5 S ( E S ( ( 5 9 1 9 0 E C S ( 5 0 9 5 ( 9 5 5 6 6 ( C ( ( 5 5 5 1 6 9 6 9 9 ( ( [ 

n mnnnnn mini 1 n i ; 1 1 n 1 1 1 1 1 m mmii nninninni'inmnni 

( M M Ml Ml M It mini II I M II 1 1 II I Cl 1 1 lie In I II I II lllllll 1 1 1 1 1 M Mill 
MM3M9MI llliminilliniiligi!!limmil!!m3i!!i;i|))!linimmH!l!!! 

I I 1 ( I t 1 I I '» » I) i l» >1 II H I| II fl it U ii H 0 fl I' fl Ii M |! M }) Jl JJ <1 ;) ;J ij <i U n <1 U |( |> || || ,, ,j •; ,> jj u jj j, j; j, |, l( ( > t j | ,, „ ,, ,, , t (] 


DT time step (seconds) 

OT time interval at which output occurs (seconds) 

ET length of simulation (seconds) 

TIM time interval from t = 0 at which the colder or warmer surface 
temperature reaches its desired value (seconds) 

RF fraction of surface black body radiation used as net upward 
infrared flux through upper boundary 
KW a constant multiplied times liquid water content raised to the 
two-thirds power 

TP temperature in degrees Kelvin used to initialize the temperature 
profile to isothermal 
ZO roughness length (cm) 

DTEMI change from surface potential temperature in column 1 needed 

to obtain a uniform surface temperature among the other columns 
when a uniform discontinuity in surface temperature is desired 
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UF friction velocity which enters into Liu* quasi-adiabat ic 

velocity profiles when they are used to initialize the veluci 
field (em/sec.) 


12 


b. Control Index Card 


IPT 


lit 

i*ii 


m 


i 

IRIS 


Ml 

I ii i 

I I 

in 


i 


in') 

Mil 

IS'iS 

ill! 
Dll 
1111 
i n i 

i*iii 


IW 


1 1 


I SCO 


s a o a( 3 : o b s 

I i ;l 


I I I I 


Inn in 2 j 


I RAD 


minimum 

iiiiimuiiu 

55555555555555 

iiiiiiiismii 

n n n n n i n i 

I I III 1 1 1 1 m 1 1 

l!)l!lll]|liil 

i ii n 'I ii it ii ii < ii ii :i :i it 


1RSFC 


i c|c n ii cj 
1 1 1* 

i |i 1 1 1 1 1 

nlnn ? 


IU 


one] 

•llll 

I t I 

mil 


nm.M 
con c c' 

1 1 ,i ii i 

i n i i 


i 

IP 


|3 13 0 5 

• n ii (l , 

III I 


minimi 


II 

,n c s n 

■ ii ii .* 

mi n 
in n n 


i 

KT 


I 

KR 


15 

NK 


1 


3 0 0 0 t|3 0 I Cf! 0 0 3 5 

lull 

II llll ill I 1 1 

.! 


IZ I MOD 


csoo; 

1 * '• 1 i 

1 1 1 ii 


mnjnnjnmtnm 1 
ijiiiiinmlmnmmimmmmninmii! 

I I II 1 4 4 4 4 4 4 4 4 4 1 1 4 4 4 4 4 4 4 4 4 1 1 ' 4 4 I 4 4 4 4 I 4 4 4 4 4 4 4 I 4 4 I 
5 5 5 5 5 5 J 5 i 5 5 5 5 5 5 J 5 5 5 5 5 5 5 5 5 5 S 5 5 5 1 5 5 5 5 5 5 1 1 5 1 1 5 5 5 S 
iiiiiiiinmiiimsiiimtiiiii!; uncut nice 
iimmimmmii'iimmmmmmmiii 
1 1 1 1 it n i s j e 1 1 ( 1 1 1 1 1 1 1 1 1 1 1 1 ( 1 1 1 1 1 1 1 1 1 1 $ h o 1 1 1 n 
n m 1 3 1 « n : n n ; o s s n s i n i n s s s i ; n : : s i s s j o 3 o 

in n a r n i* :t . 1,1 i> ■« •» i u r <i u y ji w » .1 is w 11 1 > »; i> « v it t i +■ 11 11 n 

ImIi |C«I 


semen 

ill 1 -I *i 

1 1 1 1 1 1 
mini 

1 9 3 3 3 1 1 
I I I I 4 4 I 
15 15 5*1 
(Of E(C( 

mini 

m m 1 

9 9 9 9 9 1 3 

1 ij 1 1 1 1 11 


IPT = -1 
= 0 
= 1 

IRR =0 
= 1 

IW =0 
= I 


ISED = 1 
^ 1 

I RAD = I 
+ 1 

IRSFC = 0 
= 1 


IU => 0 
= I 


nWT ,,, t , ,1000. .286 . , 

PT(I,K) = 1 P ( — "p — ) , isothermal 

PT(I,K) = PTI, adiabatic 

PT(I,K) = input data, based on field observation 
R(I,K) = RI 

R(I,K) = input data, based on field observation 
W(I,K) = 0 

W(I,K) = input data, based on field observation 

with sedimentation 

without sedimentation 

with radiational cooling 

without radiational cooling 

R(I,L) = R(I,2) 

R(I,l; == saturation mixing ratio at surface , temperature 
and 1000 mb 

Quasi-adiabatic U and V profiles 
U and V = input data, based on field observation 
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IDTKM 


HR 

COOL 

IP 

II 

KT 


KR 

IZ 

NK 

IMOD 


* 0 DTIiM(I) “ 0.0 

= 1 DTEM(I) = DTHMI between ITEML and ITliMR, and equal:) zero 
elsewhere 

* -1 DTEM(I) = input data, desired ground temperature distribu- 

tion 

,_2 -1 

= net radiative flux, cal-cm -min (positive upward) 

= rate of radiative cooling, °c-hr * 

= 0 no HR and COOL output provided 
^ 0 HR and COOL output provided 

= 1 turbulent exchange coefficient does not depend on local 
condition 

= 2 turbulent exchange coefficient depends on local, conditions 
= 1 two-dimensional simulation is initialized from list 
generated by one-dimensional simulation 
f 1 vertical grid level to which first potential temperature 
in list is to be assigned 

= perform the same function for mixing ratio and liquid water 
content that KT does for the potential temperature 
= specifies the number of the column for which wind, potential 
temperature, mixing ratio and liquid water will be punched out 
= vertical grid level to which hygroscopic material is to be 
spi .i 

= 0 . odifieation by heating 

= 1 modification by spraying hygroscopic material 
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c. Grid Specification Card 


* 1 • 00 ftp ♦ 1 

o.2ooe + o 

1 . 00 1 C « 0 

1 0 

30 

1 

e 

1 

itoh. 

0 

ZAI, 

ZAK 

XAI 

IE 

KE 

II. 

IR 

IThMR 

1111 111 

1 11 91 

9 9 9 0 9 

1 1 1 

9 0 9 i 

tote 

uc; 

HOC 

T p 1 2 

i i l i i i i i i 1 1 

il i| ii ,i ii ii ii i •): 

a ii ii i n !■ :i a 'i 

nan a 

a a a i i 

It 1! Il 1) l 

i' 


. . 51 , 

11 1 1 1 1 1 1 

111 1 1 1 1 1 1 

1 1 1 1 i 1 1 1 

1 1 1 

mi 

1 1 1 

1 1 1 1 

1 1 1 

Ml | 

1333333321 

232 32331 

322332222 

2 2 2 2 

2 2 2 2 

2 2 2 2; 

2 2 2 2 

2 2 2 2 

2 2 2 2 2 


] i 
i < 

s s 
u 
1 1 
1 1 
t 3 
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3333333 133333313 113333 
444444444444444444444444 
S 353333333 535535555 55 
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lllltil ilMIMU imu 
3 3 M ! J 53 1 1 ! J i S 3 U 5 5 5 ! ! S 1 

t I I II 'I II 1} (I It it II II II II II II ,i II II II II II II H 


1 0 0 0 . 


nrxx 

3 tcoceccoocsooca 

II i: II I. II III’ '1 'I i II ’I I i i i |i 

1 I I I I I II I I ! I 11 I I I I 
2222222222222222222 
3313333 13313333333133333131111311 111311333311]] 
4444444444444444444444444444444444444444444444444 
5 5 5 5 5 5 5 3 5 5 5 5 5 5 5 5 3 5 5 5 5 3 5 5 ‘ f 5 5 5 5 5 ; 5 S 3 3 5 5 5 5 3 5 5 3 5 3 3 5 5 

E S S S S S S 5 E E I E i E £ EE 5 E E E 6 S E S £ E E E E E E S t ( ( E E G E E C E t B £ £ l E 

iiniimiiimmiimiimmnmmiiiiiiim 

I I I I 1 11 1 4 t 1 11 it t 1 1 t I I I I I t 5 I E t E 3 3 1 I 1 1 1 1 1 1 1 M 3 I C 8 

I IS 13 ! ! 5 313 1! 3 3 ! 3 H 5 3 1 ! 3 ! ! ) 3 3 S S 0 3 5 5 1 1 U 8 S 3 9 3 5 5 3! ! 

> II II H II 11 1) II It II ii il il M ii il li il il l) li II s' :i ii 4 i .1 :] l, |i li il ;i n l. |i if ll I 1.1 II . i| ■ i n 
• Ml* ICO 


ZAL = height of lowest vertical grid point in the atmosphere (cm) 

ZAK = expansion factor in the vertical 

XAI = expansion factor in the horizontal 

IE = number of horizontal columns 

KE = number of vertical levels 

XL = leftmost grid point of unexpanded grid 

IR = rightmost grid point of unexpanded grid 

ITEML = left I-limit on non-zero, uniform DTEM(I) 

ITEMR = right I-limit on non-zero, uniform DTEM(I) 

DELX = uniform X-spacing (meters) 
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Output Formats 


UH tFCB-2* lOtU.M 

kick i i c s 1 1 ! o o n c o e s e c s : a o c c : c a 5 s o 0 a c 9 a o o o c c a v a o o a 3 c c a a a a a r c c c o a a * a c c . • a a ; a 

I I )•« I f ’ M <1 " 'i <1 ’l 1 l l l ;j ' l: u 1 1' l . I :• . i .» * I i a ‘t # »' *' *i M «« i u »» .t: . .* U . • I » »• i • »' * • • • *' 1 * » 

I It 11 1 1 11 1 1 1 1 n i II II 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l i 1 1 i n 1 1 I i ■ i 1 1 1 1 1 1 1 : i ; i 1 1 n i M i i 
limmumii i n i i i ii!iii)iiin!m!ninmmmini!i!ii!m!a:!:i , iiii!: 
iiinmiiiiii s jniiJ ijiiijiinnnni»)ij:>nmijijjij::n 3 i»j 3 :jjjj: ana 
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5 5 S 5 5 8 i i i 5 5 5 5 5 5 S 5 5 5 S 5 5 5 5 5 5 j 5 5 S 5 S 5 5 5 5 5 5 i 3 5 S i 5 5 5 5 5 i 5 5 5 S 1 5 5 5 . S 5 5 5 5 a 5 8 G 9 5 5 5 1 ‘ 5 5 8 5 
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m m m m a m m 8 m i m m us s m s 9 5 m 3 s n a ! s a a m s j » : im a m t » a m a a a ? - « j ' . 

I I I i I » 1 i » u i> J i* « u tl :mj u ii b :r: :i ,*:i r ii ,i j n uuiu* ji 11 »< i> •’ mm •' « » *t « r »: •* » * i » ». i r . .< ». t.i :# •• i >4 

(IH Iptae.M OPlDMZx.) FMTT 

0 8 1 0 0 0 0 1 0 0 to C * - 0 0 3 3 C 5 0 C 0 0 0 0 0 0 0 0 0 0 £ C 8 » 0 G 0 0 3 0 S S 0 3 9 C C 0 C C C C 3 0 1 0 3 C 0 0 C 0 1 J C C C 0 J,; 

1 I >*ll 1 I l >1 imi oi i n •! I .»;i :: tt ,» a u i| :i m it :: jt u r- n i> :i *1 n r *' »j u « <* « <i • u • :: u * » i > ;i 'i ») 1. 1? iu< »? n i ll t. i , t , t • u 
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_ l*.u«»C»i 

The printed output for the various quantities is in convenient 
standard units as shown in the format heading statements 3000-3800 in 
the program listing in Appendix C. 

FMTH = Format for hieght and variable informations, including wind 
velocity, R, liquid water content, RAD and visibility. 

FMTT = Format for height and temperature information. 


e. Microphysics Parameters Control Card 
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11111)1111113 111) 1 1 I 1 1 1 1 1 1 1 1 1 1 1 1 111 1 1) 1 1 1 1 1 1 1 I I 1 1 1 II I I 1 I I I 
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<Mi) ict> 


NH = number density of hygroscopic material for modification 

VHF = Van't Hoff factor 

PDEN = density of pollutant 

HMW = molecular weight of pollutant 

VHFA = Van't Hoff factor of the other pollutant A 

PDENA = density of the other pollutant A 

HMWA = molecular weight of the other pollutant A 

RAHI = initial radius size of hygroscopic material 

HDE^ = density of hygroscopic material 

VHFM = Van't Hoff factor of hygroscopic material tor modification 
HMWM = molecular weight of hygroscopic material for modification 
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f. Aerosol Particle Control Cards 


NDIS(l) N 1 )IS( 2 / .*MIi>( 3 } 

lEEICMi' 0 C 3 fl » 0 S Mtjjltl C C « 0 ICC J 3 8 9 t 6 E 5 * t C C C 0 t C » 3 C C C I 9 5 3 6 C 8 0 I 9 CC CC J I C , ,1 BE 
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4<<<MIIUHM4< M4M<MM44MMM4<4M4<M4MMM4MM<MMM<MMMMMM<M4 
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g. Other Input Parameters 
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RI = initial value of water vapor mixing ratio, if .set IRR = 0 

PTI = initial value of potential temperature , if set IPX = 0 

OPTIM = star*; time of modification 

HEAT = source temperature, if set IMOD = 0 

YI = initial vertical wind velocity 

DTF = time step factor 

OTF = output time factor 
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li Variable input, Non-uniform wita lloihpl 


1. Horizontal Wind Velocity (X-Direction) 
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2. Horizontal Wind Velocity (Y-Uirection) 
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3. Temperature 
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4 , Water Vapor Mixing Ratio 
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Liquid Water Mixing Ratio 
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B.'s. Lint of Computer Variable . 

The following list of variable:; is not arranged in alphabet, ieai 
order, but rather variable:; of a common type of variable:; which art 
concerned with a Riven prop. ram operation are listed fogethei . 

a. Grid Specification 

1 . Vert ical ( h’id 

K - index for vertical grid 

KK -■ number of vertical levels 

KN = KB - 1 

ZA(K) - vertical coordinate at level K 

DZA(K) = grid spacing between level K and K - 1 

ZAL - height of lowest vertical grid point in the atmosphere 

ZAK = expansion factor in the vertical 

2 . Ho rizontal Grid 

I ~ index for horizontal grid 
IK 3 number of horizontal columns 
IN 3 IE - 1 

X(I) ~ Horizontal, coordinate at column I 
UX{I) - grid spacing between columns I and I - 1 
DKLX ~ uniform X-spacing 

IL 3 leftmost grid point of unexpanded grid 
IR 3 rightmost grid point of unexpanded grid 

XAI = expansion factor in the horizontal 

b. Variables 

1 • Prognostic 

PT(I,K) = potential temperature 
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RW « gas constanc for water vapor 
SIv’iMA ■ Stefan-Boltzmann constant 
2. Variable 

KW * mean mass absorption coefficient of fog for infrared radiation 
(cm 2 g -1 ) 

RF = fraction of surface black body radiation used as net upward 
infrared flux through upper boundary 
UF » friction velocity used for computing KA 
ZO = roughness length 
CV = sedimentation constant 
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15.5. Sample;} of Output 

1. Temerature 

2. Relative Humidity 

3. Visibility 


Three samples of each paramerer are illustrated in ulie next t lire 
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IMf tiff A T * l?n o i 0 (jl) e In 1 

I! Of >■ TiCf. uf.L TIM CTfP Af«i* COVPuTt Nf » pr v . ;>,Oi T ; C : 
V M> IAM'.O 

SU. F>l>!iTJ*. - tT . » 


F A T' 

f A L * <i ■ 

D I m '« , ; 
CO" mo,/. 


CO ’“••O’, / 
CO" • Of;/ 
Cr i •„/ 
c r ” •* f ’. / 

C o m f c f, / < 

c<> •' <:>, / 
'/*"(••,/• 
C 0 ! I v 0 f, / i’ 
COMO’./ 

k 

c o M " o *, / 

k 

CO 

if err- ; . 
Ifni",, 
i f c t r • . 

I F CJK , L 
f, t> T - *j 


n „ > - 1 v i 0 7 - ’ 1 , J V s 1 
k* ( I f f I u I ’ 

7FnCj/JiFOC(jX,j2) 

■ i T F C ( J / , J i , J r > , TfOACjn^/^r) 

> l, I, / » / ( J a , J I i ,1k ( J < , J Z ) , J N T < j * , J ; ) , r » T C j 

t>kH(jX,j;),VIflTY(JX|J4> l %Ulj(JY)|lkC 
,POl vAC J < ) , fJSF , 

0 1/ F U ( J 7 ) , F u ( J 7 i 

u'/ o c j * , j ; ) , c v ( j ;• ) , f v < j : ) 

1 i / v c j a , j ; ) , r r c j : > , fy(j?j 

o a / i< c j * , j : ) , r i> < j > , f » < .j ; ) 

t. 5 / u ( j * i j , f » c j ; ) , 

V. * / F U > a , j ; ) , t . t>T ( J 7 > , F f T ( j 7 > , T < „ X , J .; ) 

J - 1 n j * , j / ) , (j;) T‘ :'i jx 5 , . x < jx) , *• < j r 

0 M <.,*(Jfc)| r (Jx)kIt.tf-f'kTI"i, ■ it » i 2 F L ’0 i 

n '» / 0 , c f’ , 0 , 1' T , T I v . , C V , 7 A f- , C l , C •• , C 1 , C f , F , * 
1 1, 0 f C , I I . ' P T I " k H 4 A T , 7 0 , C tk , £' F N , P • T 

Lin/ FCloX|Jt)|*Aff(JX|J7)|KA|, Ijf,i i N M i H 1 1 

■1 M f A , ( ■ 'j r \ A , p r X A , H f) l f; , V I, » !' , r, " „ " 

oil/ • AC C Jx ,.! 7 , JY> (JX , J7 , JY) ,FV<J 7) 

o f. . :• t ) c . t -• [■ r 
. o. -n a . > p t i y j r t j 0 <- £' t 
0 ! . < f T ; J .,0 Ty 1 J „ 

.17 U; TJ l.’O 


X i J 7 5 • 
i 0 t J Y ) 


tif t J 

>MI 


* i J 7 ) 

0 A ( J r ) 


5 » o 7 C J 
i T (. " » , 
•< . I M I 


I r f , 

0 > U , 


7 a Hi) 
MF 
I RAO , 


,hT , Vhf , F t : , H ^ , 


I i r - o 

: r ( IV it, .. I, . , ) Tv 173 
I ) C 1 2 0 I '. , ' 

7-0 120 <n:,,i,-1 

X F ( T V't kt . / yPT I ) F «»1 ( I , *.K ) = K AMI 

IF ( X . f. , f. f ) go T ( ;t 

I F C k , 0 i . . 1 ) r,g TO 7y 
F C C ( I , > » 1 ) - . 

K A M ( I , K * 1 ) - f> . 

F C C ( 1 , K ) = 0 . 

IkPAHsPAHlI ,F ) 

H AF’r A * 7 , T A 1 / * F, A li i * * 7 * ft I *" f. / 5 • 

CALL B ». f C F A H ( l , f ) , I (I , X ) , P c I , K ) , P ( I , ». ) , V H F ", , ,i f; ' , t < •• „ ' , f.h ) 

TPAMsF.AH C I , X) 

T V - 1 . 2 0 I * ! * < A H C 2 , y ) * » 7 
0 1C T = 7 A {;,<) -i A ( ’.X-1 ) 

F C C ( I , L ) - A* 7 , 1 - 1 6 * CXAMC I , F ) * * :-PF)An- * ' ) *I,H/ (I *o ' L T ) 
ft AM ( 1 , k. -1 ) = TnAl! 

COL TJL.IF 
CONTI?,, II 


9 B 


I J OOK 


0 T r !• I !, I / I \ 

f'PT-Tlr MCT-l't T 

If CfT.Lr.DDT) u U TO 130 

NPT S CT/M.T 

DT = DDT 

FCCCI |K)“fLCCI t K>/UDT 
130 CONTINUE 
C 

C UPWIND r ’D Lf’UP OVTu HORIZONTAL ARID 3 Y 3 T l V. 

C 

DO . 0 I - t , i * 

DDX?DT/D*(l) 

C 

C DOWNWARD Dt, LOOT 0 , f i< VIRTIChL On I D 3YSTIK TU 3 " T UP IMPLICIT 

L INTf oRAHO.. 

I 

DP ?0 N--0,n’. 

K-KN+O-N 

l 

t CONFUTE wOnHNu VARIABLES 

C 

JFCTIR'i.LU ,DT) TPR (K ) sPR (K ) 

C D = D T / 0 Z O. ) 

D P Z = C P ♦ D .. A ( c ) 

Arft «KA < 1 +1 ) / < 1 . + ZAK> 

AR = A 

Alls ! D*K"C 1 , K + 1 ) / < 1 . *Z AK ) 

MIU - : l * A u 
A T = A 

C = DD*KACi,r. ) 

CW = C 

C IIS I'D ‘ C 1 ,K) 

CUU-II*Cd 
C T = C 

IF ClSfO ,Ki. .1 ) 00 TO 10 
l 

t DROP SlDIA’cNTAfluN 

C 

AW = AnK.»*«LVn( h (X,, <+ 1)*».A6 7)/L. 

CV = L n-l>t Z *L V« ( n (I , K-1 ) * * . 067) / 0. 

C 

n u = i . * a ♦ f 

uUl<= 1 . ♦ AdL t cuu 
L T s 1 , * A T . C r 

D P T - P T (1 t K ) -U ( X >0'* (PTCI ,0-FTCI-l , K ) )*PuX*C1 .1 , n . / P C 1 , K ) ) **,^ a s* 

1 L * f C Cl ,n )*M/CP-Y(I ,K) * CPTCI ,K) -OT(I , N-1 ) ) HT/PZA ( Y ) 

I F C ( l v CP. U . 1 ) . AND . CTIMF .GT .OPT t") ,A(.D. (K.LO .NO > DPT = DtT + (100i)./ 

* P(I,K)).*J.Dio*L*FCCCI,K)*DT/CP 

L 

C RADI ATI ONAl COOLINfa 

C 

1FCCIRAP . LO . 1) .AND. C W ( I , K. > .CT. U.J>> L,PT=DP1“CP*W(I,K)*CPTCI, 
11)**A.)«iXF(-I,<T(I,kL) + INTCl,L))*PRCN)*0T/CW(I,Y)«*,.'T3> 

C 

DU = U C I , K) -L C I , r.) * CU Ci , K )-U C 1-1 , K) ) *AD»- (F C i , R ) - r C I -1 , K ) ) - ODX/DtN+ 

1 F * V C I , R ) « D f - Y C I , k. ) » C U C I , K ) -II C I , -1 ) ) * PT / 0 Z A C k ) , 

DV = VC I ,K j-uCl ,r.)* <V Ci ,X.)-VCI-1,IO)*PUXvF»IuCl,t.r)-uCl,k,)>*DT 


r \r.E ® 

quality 


1 -y(i, c, o-vU|i'-i)>*ri/&i«( ►. > 

DW~>. ( I , « )-t.Ci ,M * U ci ,► 5-WU-1 , O ) »f I, * + f C (I ,K ) ."T 

i - v ( i , «, ) k c » c i , > ) - a i , y - 1 ) ) u r / c : a < ► ) 
t t(=r, t : , * > -i, a > . u t; ,k >-p < i-i , ► > > . I x~i c i :,*■>* k t 
\ -t u , o • c < a ,o-i c i ,k-i ) > •pt/ccaf* > 

I F { ( ;•• ,.c . . 1 ) .AGO. (TI'-T .0 T .OFT IM> .Ai.C. Cr .N*. > JOR-GR-I CC C I i K> *l‘T 

t W{»)‘Ck/(i -*»•{»(*. *15) 

f.llCr, } s ;nu/ t'll-A Jl!«f II O' *1 > 
i VI* ) ■ 0-G('-A„fl» r V (K *1 / ) 

tPT(x)=,f/lT-AT»t|T(»+ 1 )) 
i <' ( * ) - C / ( 1 - A • t h C r. + 1 > ) 

mi <* ; s ci, + «';u* mu ♦ i ) ) » mu ) /ryu 
f v ( * ) « i v + # ru » f v i 1 > ♦ i ) ) * r v ( *. ) / f : uu 
f l> T ( K ) r I * ,,T * f F T l K ♦ 1 ) ) * r FT (K ) / CT 
F P ( K > ■ : ( p I + * « f r.t >: ♦ 1 ) ) » f ,J { K ) / C 
P f (X ) - k"1 1 , K+ 1 ) * ( Y( I ,K*1 )-Y (I ,1 ) ) 

T T " T ( 1 , 1 J / I F K ( K ) 

P C I , K) " p '. I,k-1)-(lJ(I,i')*(Y(I t K)-y(.-1 I K))/l)''U)*YU|K)*(Y(I,K5- 

1 -PP ( K ) / ( C i ( K ) « ( 1 •• l A., ) 

2 > n " C ,s )* (v u ,H)-Y (I , K-1 ))/[•; O.) *C* ( Tll.Fl/ TT -1.J5* 

7 t> t. f. ‘ :> i a c ►. ) 

2u COM IN CM 

c 

C UPDATE So SPACE ir i UN t> « ft Y CONDITION 

C 

I F CTI*t at . TIP) PTC I, 1 ) =PT<1 , 1 ) ♦CTC'! J) .TICf/Tl" 

i f c a •■•<>;> . r >. ,c) . a.,p. n i t- t . i- t. op Tin .an;?. c i .or . a> . and . < i . it . M ) 

* FTU.D-T'.M 

t c i , i ) -n c i ,i ) 

IfCIR.rC ... ’i. J) ”( 1 , 1 ) = FPC2) / (1 ,-Lli < n ) 
l 

c 

C SATURATION '.or. PACT PoUNPARY C 01, C I T I 01 . UN P 

C 

I F ( I K 5 I C • l « • 1 ) l: ( I , 1 ) = P S F ( T Cl , 1 ) , 1 0 C . ) 

C 

C Of>»Afcp U LUl’ 0,PF VPPTICAL CSIt S Y L T £. « T TL COKP'UTfc Nf * PROGNOSTIC 

l AND 0 1 A i.OGTIC VARIA.iL f $ 

C 

I N T C I , 1 ) - 0 . 0 
CO 7 0 F = s , i t 

ll(I 1 k)=u.(i.)’ C 1 , K - 1 ) * II C > ) 

V ( 1 , K ) = v C N > » V ( 1 . < - 1 > ♦ F V ( >' ) 

rII,A)'TU,i.-))-UKI|ii)-utM 1 i())tjHU)/cx(:) 

PTCl,A)=iFT(K)*Fr(I,r. -1)*rPTOO 
R C I , K) = ' I. C A ) * f ( I , K - 1 ) ♦ F R ( K ) 

PR <* ) = (100. ./FC I ,K) ) i 

T II ,K) sr-T (] ,k ) /FUCiC) 

T W - W ( I , K ) 

k(I,ll)=U(K)*,(l,l,-1)FlklF) 

t 

c computl i. a t i u o of n a r i ’< droplit nr uiuii «atfr co.ntcnt 
c 

IFC*(I,K).LT. C.O)*<I.K)»P.O 
D M ASE=(^(I,I. )-Tr)/A<-u{FC(I,k)*OT) 

DO 22 J -1 . j r 
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I 


I 


i 


t! 

*| 

f 


i^s 


A w A r - C = t> ' / S u * I F C ( 1 , f. , J ) 

1 f C.-5P.' t. . 1 ) o Tf' 23 
A’"ACA = S v a:^ »TFC A ( 1 i K ) J ) 

T? bTAsA 'i.Sa' 3/ ( A* . 7 . HU*HGI«-A( J ) > 4 ^ A £1 <1 ,fc , J ) . . J 

TI?BI5A = PAE,(IiK.(J)**? 

IF(TESTA.LT.TF.blSA) TF.GTArrTIUISA 
R A f A • , k- j)?TISTA-*(1./3.) 

<.3 toil T • : ■ 1 ■ 

T E t> * ■ • ?/ (A*i. 1 41 6*NDI S ( J , ) +hA0(i , K , J ) * *3 

TRIM' - ' I i K i J ) * * 3 

I F < T ' , . T [t f> X 5 ) TCbT = TRDIS 

H2 ftA:;i,<,j)=TLST»»n./j.) 

t 

C COMPUTE hF.y, VALOtS OF 1 NT , Co T i AM) a 

C 

I M T ( I , A ) s I f, T ( I , - - 1 ) ♦ C 1 * ( U ( I , K ) * * , (< 6 7 * * { I , K - 1 ) * * . c *■ l ) * 0 1 A < K) 
l,w: = E'l.iT((l.U,K)-U(I,K-1))**L , f(v(I,K)-V(J,K-l>)**2+(Y(i,K)-Y(l,K-1 
*))**2)/D..AO) 

I F ( W ( I , K ) . CT . 0.0) 60 TO 5 0 
CPT(I,K>=CF 

t 

C D I A f N 0 5 I i OF f ACHAf.GL COEFFICIENTS COVlRIMG M EllT k A t . i T A r L E AN D 

C UNSTAPL" C 0 M> I T I 0 N 3 ijSTWEF.« K = 2 AIIU THfc bURFACL 

C 

S=(PT(l,n)-n (I ,K-1 )) /OZAU) 
i f < k . r .. 2) ib = s 
GO TO 60 

30 CPT(I,C)=CP+tH*HbF<T(I,lO,P(I,».))/(T(I,K)*»2) 

5=(T(I,K)-T(I,it-1))/0ZA(K> + CS/CPT<l,K) 

I F O, «t()» t) S b - 6 
1FO .11, L ) 60 TO 71 

i f f s s .or. r.) go to 75 
I f ( !< S . Fu. r.) V- TO 77 
PHt = liS.RT(or.-RT(1.l/0-1b.O0*2(K)/TL)) 

K» (I , K) s . « ' l (I. ) *U» 1 » PHI 
K A C I , K ) = a ( 1 , K ) * P H I 
00 TO 7 -7 

73 Kr(I,K)=C...*tU)«Uf1W(1. + 5.»U<iO/TL>> 

GO TO 7’ 

77 K v ( I |K)=.A*t(K)*UF1 
7 0 y. A< J ,K) -aM( I , A) 

C-0 TO 7 7 

71 FM(I l 2)*.H*»i.KT(UU,2)**?*V<I,Z)**2*Y(I,2)*»’)«tZA(i)/<AlOOUA(2> 

♦ / 7 0 ) ♦ * 2 ) 

F. A ( I , 2 ) = K M (1 , u ) 

U F 1 s S 0 R T ( K ( I , 2 ) * p U Z ) 

I F < r .to. C.O) 60 TO 70 
TL-(UF1<*3)*PT(I,2)/(iA*G i KA(I,2)*!)) 

GO TO 70 
i 

L SHIR-TYPl l XkOMNTI At UGCRIASF 

C 

79 DfC=lXP(-?.'F*i(A)/UF1) 

AF'<I,K.)=AH(l,K)*O r C 
K A C I , K ) =AA< I , A ) *DFC 
L 

70 CONTINUE 


iff 


OR 1 ’- ' 

QF FOUR QUALiIm 
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L 

C OP OAT & '■.IMU •. «•.**»» APT C«’JMTION O'; a Af HP. 5 A T ul' A T I C l» 1 0 J " - T ;•• t N T 

0 IF NO FLUX CO'aix I T It; N <„FD 

C 

if < x n g f c . l . o ) n < i , i ) = r t » , ) 

0 

PO C f N T I Mi f 

0 

L COMUTt .r.vXTH OF »» A T :: F. OPUPLfT 1 Y w i f AC I’ll Y $ I C 0 ..OuATIO' 

1 

0 0 ! I - 1 , 1 

DO JS X " i x f 

F^ISftlU.Fl.FCiK)) 

AH< I , X X "it (i ,F. ) /Pb 
TVl'.-O. 
f C < I , A ) * w . 

CO 74 J = 1 , j Y 
hPAC = SADCI|X iJ) 

U5f = 4 • ? . 1ft U«ltl IS ( J) **■»*!>[, -F,/ 5 . 

call 1.1 F U»t ( I , K , j ) , T C I ,k. ) ,F'C1 , x) , a (I ,x> ,VHf , 0 I Z" , f< C 1 5 ( J X X 

TFC (I ,F , „ )-ft‘ ' . nit • (PA DCI , F , J ) “7~r i,AC.»*2) *f| US CJ X a. 

F C U I ► X = F r ( I , X ) ♦ I F t. ( 1 , r , J ) / c ; f 0 * c T ) 

TVIS = TVIO+(.CIbU)*(IAOCI,K,J)**’) 
lF(Mr,'„.1) <0 TO li 

H P A r e T A r ( I , F , J ) 

0 I = ft * 7 . 1 ft 1 1 ’ i< i. i i .«( J )**.’■* P C ; A / I . 

call t. X. F i r< » i: 1 1 , F , J ) , T ( I , k. ) , p ( I , F ) , l( C 1 , F > , v M f A , C I S f‘ , m a , M’ I r . « ( J ) ) 

TFCA( I , t , J) ~U-: . HI '• Cl. At- ( I ,F , J X ** '-«AAF «*i) - fiU I OA (J ) / i . 

FC(l,K)=FCU,F)nFLA(I,K,J)/(Cr.N*r>T> 

T V I S ~ T V I 0 + !, D 1 0 A ( J X * ( H A t ( I , X , J X « * 2 X 
j 7 C 0 NT I *. LI • 

54 CO’lTI’lUC 

TVIS = Tl/Ib + M'«CrtAnCl, K )** t >) 

V I b IT Y ( I , X ) ' i . -r 1 2 / < ?. ‘ 7 . 1 4 1 6 * T V I E * 1 J C . X 
5 7 C P N T I N U ii 

1 \ t = i •, r r + 1 

I P c i ■< r T > L C < A C T } T I !• f i T I v f ♦ o T 

I F C : M> T . i. 1 . M, T X ,0 TO 12 0 

1 F ( V. . L 1 . 1 ) 0 T - L 0 T 
1 F (TI’U >V ,T . CF’T IlO f,K = NK-1 
A T l.i P t. 
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original 

OF POOR 


4 

l 

t f'l'lf.T j::. l-vi rif.t 

c 

SU'.'PCUT I . C I h f. T (.!, f Oi,M) 

C 

PABA.Vi.Tf i, jx-l J , j : = » 1 , j Y = 1 
envrvK/i L?/ Z6Ui),XUX),It.,Kr,TlM ( .,DTiV|, 1T( .v L ,i Tt fi 
0I"f NSIC'-. U(JX, Jt) ,fuR.'(11) 

L 

C OUTPUT D 4i C 4. !' L M A T 1 4 , 1, 

C 

a E " x f - 1 

nl ; ITLL f i.i;/j) > I ■' f I r T £ V I , J T p M L , I T £ f R 

V. 

C PA'jT. 1 Cul. i, PA 3 1-10 

C 

(XII), 1 = 1, 10) 

on 430 J-1 ,* L 
K~r r -j +1 

430 W»ITLC6,fOi, •■■*) OA(K) , <0< I ,K) ,1 = 1 , 10) 

L 


I F ( I t 

.LT. 

1 1 ) 

0 u 1C 4 V 0 

PA. 

4 C 41 L 41.' !. 3 

1 1 ~c 0 

«r;tl 

( £ , 4 1 V 

S) 

(X ( I ) , 1=11 , 00) 

DC 4 i 

C J = 1, 

r K 1. 


K = r. C - 

J ♦ 1 




460 .RITl (6, fCht') 4AU),(0(I,K>, 1 = 1 1 , 2 u ) 

4 

IF(It .LT. CD LO TO 4 VO 

4 

C PACF 3 CuLUHvS cl-JO 

c 

» p I T [ ( i, m 1 0 f ) (x(I) ,1=01 ,30) 

DO 6?r J = 1 , I L 
>. = K ! - J ♦ 1 

4?U kP IT t ( t , mjk: ) 3 A t K ) , f 0 C I , F ) , I = f 1 , 30 ) 

4 

I F ( It . L T . 31) UU TO 4 90 
4 

C PAC-t 4 C 4 L o '■ N S 31-11 

4 

* R I T 1 ( 4 , 4 1 4 3 ) < * ( I ) , I = 3 1 , I 4 ) 

DO 4 4 0 J = 1 , r. L 
K = KF -J + 1 

400 WRITK&.FCM*) 3A (K) , <0(1 ,K> ,1=31 , U> 

C 

490 COMIUUO 
KE = LL ♦ 1 

4103 FOP"AT(1li1,(hy<Yv 1 > a ,S*,10(- r PFl?.?) / 1 M , i> , Sli? It'D ) 
400^ rouKA T ('1.01 I (■• .. - f f V , 2 , 4»l SI C , c X ,',hpti f; I = , f 5. 1 , ,’n i, 4 x, 
lx, A»< jt i: , u) 

4100 f CI.MAT ( 1 H 0 , (HA U K) ~ , 10 (-3PF 1 D . 0 > / 111 , 3 X , c H .' ( f ," ) ) 

R r TURN 1 

LAD 


1>A.GE IS 

QUALIFY 


X II I T F ML 
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¥ 


1 


c 

i 

i 

c 

l 

L 

c 

i. 

t 


l 

c 

L 


L 

l 


*<■ 


50 


t <••uf.it :, : i ■ u v »-t P„Ti'.f 
G 0 t V t i > ‘ i » T >i n n ! „ t . u '• T 1 1 N 

*• * r !• / 1 t - . ■ ( .-«/!•♦ / . * < 7 5 
GU.'I C t' r I *. t (' f. c , T , t', it, v»it ,»> >' 


►'h y j c ,.l c u *■ g r .»*. i g 

T ■■ 

Ai.f HA".; , w ? 

i f r a - r ,■ •>. i 

v, A * A -1 . I. 

A V ' ° i 1 1 

M-r/./ „ 

C A Y ‘ t. T f a i Y 4 % ij 4 C ■ » G 7 A 7 G 

r- r : ■ g k ; , i j 

l - 1 1 . i 1 1 * t * * i . / *. , «• 

T 1 • .'.,.*1 

t r ' ' . 1 : T - 1 'JO. G ’-’Gi - A 

: i ( - ’ j . ;■ - . . i .. * v t - r ; 1 . i > 

Go* •; •• . ' •< ’i -a * i. / 1 . u / 1 ) * *r . . r ’ : n fn<t <m,. ;• n , r .• ’ * 1 1 . / 1 r » . n-i . , : n 

. ' • 1‘ *, * l' « U .a ?-;/T-G.;’’7?K<.)/r 

I. f f i - 1 . / > * ' * n. / c r. 

0 - • i ■ • ./it ii 

At 11- li, I / w . > * t ►■ / tAui’MA»IO * (;A"A~1 1 1 , 1 *. f /,, ,*)•», /( 

■' 1 , ) 

1 L - i 1 - I A / , )*(;/; , T ,\ 5 * ( : * « . u 1 Y / C H V * T ) ) • * 1 . . 

Cl A (! ■!./;•♦ At. <i . • u / (.»)*! f F F 

CA-. »'jI' A/l'iV*T«».lAAt7) 

C ;5 + A . 1 * v ii f * H’, > / m » 

1NT! i '< A 1 1 C „ 

1 =■ 0 

v T . 3 1 

t r- . 

Cl". TIM, I 
K h - >’ / s ■; 

Gii r r - p »i - I . 

1 = 1 + 1 
1 P A r f; A f< 

T I +• f — T I i- u + L 1 

f (? = '.£■» I r ..jf v - C A/R/. t'+fi. / ( BAT * * 5) ) / (RAG + f L ) 

f; A p - P A ii + L 7 * f it / a A li 

t P = A *! . HI c ' CRAo “i-Tr. ao* + 5) **./ ( 3*07 r« ) 

f. = k - p r 

liRA!>*A Miti.l -7 ii A . ) 

If (I.r.T.tCjC) ol. TO 50 

I F (PI- At. ,l,T. 1 . Gl-A) 00 TO T 0 

CONTIMJF 

[, = T I! 

1. 1 T U R f J 
E.NP 


1 


i 


1 
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